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Abstract 

We prove the double bubble conjecture in the three-sphere S'^ and hyperbolic 
three-space in the cases where we can apply Hutchings theory: 
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• in S^, each enclosed volume and the complement occupy at least 10% of 
the volume of S^; 

• in H'^, the smaller volume is at least 85% that of the larger. 

A balancing argument and asymptotic analysis reduce the problem in and H'^ 
to some computer checking. The computer analysis has been designed and fully 
implemented for both spaces. 

List of Figures 

1 Plot of the Hutchings function F{v, w) where it is positive in S^. 

2 Plot of the Hutchings function F{v, w) where it is positive in H^. 

3 Some important geometrical features of the generating curve for a standard 
double bubble. 

4 The steps to show that g > h if both are concave increasing functions of 
one variable. ^ 

5 The steps needed to show that g > h on a rectangular region. ^ 

6 The triangular domain of the {v, w) where numerical analysis is 
implemented. [H 

7 The steps needed to show that g > h oil a triangular domain with 
hypotenuse along the line ui = u. |27 

8 Showing that the Hutchings function is positive on the rectangles covering 
the shaded region shows that the function is positive on the shaded region. |27 

9 An illustration of the steps taken to cover the whole domain of volumes on 



31 



10 The basic steps of the proof that F is positive for v > .85ui and w > .85v 
in H^. [32 



1. Introduction 

1.1. The double bubble conjecture in S^, H^. In March of 2002, Hutchings, 
Morgan, Ritore, and Ros [HMRR] proved that the area-minimizing way to enclose 
and separate two given volumes in is by a standard double bubble, defined as 
three spherical caps meeting in threes at 120 degree angles. In 2003, Cotton and 
Freeman extended these methods to and H^, proving that the standard double 
bubble is the most efficient way to enclose and separate two equal volumes in these 
historically important non- Euclidean spaces, with the added condition that in 
the exterior of the double bubble takes up at least 10% of the volume of S^. More 
recently, Corneli et al. [C3] have proved the double bubble conjecture in spheres 
of all dimensions provided the double bubble partitions the sphere into the nearly 
equal volumes. 

Our Main Theorem in (Theorem 16. 3p improves upon Cotton and Freeman's 
results when v, w are sufficiently large, showing that that the standard double bub- 
ble is the least-area way to enclose and separate two regions of prescribed unequal 
volumes, whenever each region and the exterior contains at least ten percent of the 
total volume. In H^, our Main Theorem l6.4| extends Cotton and Freeman's results. 
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Figure 1. Plot of the Hutchings function F{v,w) where it is 
positive in S^. 

showing double bubbles are standard whenever the smallest region has volume at 
least 0.85 that of the the larger region. Computer plots [l] and [2] show that this 
is approximately the largest range on which current methods of proof using the 
positivity of the Hutchings function can work. Indeed, the main focus of this paper 
is rigorously proving that the Hutchings function is positive on most of the regions 
where it appears positive in these plots. The plots suggest that F{v, w) is always 
positive when v is the volume of the larger region. This fact is shown in Remark 
I2.16| answering the open question of Cotton and Freeman [CF, Open Question 1.1]. 
The proof. The main difficulty is showing that both regions of a minimizing double 
bubble and the complement all have at most one component each. Given such con- 
nectedness, Cotton and Freeman [CFj extend the instability argument of [HMRR] 
to show that a minimizing double bubble must be standard. By work of Hutchings 
[H] , the first region in a minimizer of volumes v,ui is known to be connected if a 
certain "Hutchings function" F{v,w), determined implicitly by integral formulas 
for volumes of spheres and standard double bubbles (see Definition l2.2p . is positive. 
In their paper. Cotton and Freeman use asymptotic analysis and intensive compu- 
tation to verify this result along most of the line = w in and all of the line 
V = w in H^, thus proving that the standard double bubble is standard for equal 
volumes. For the two-dimensional domain of volumes {v, w) in our unequal volume 
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conjecture, the proof that F{v, u;) > is considerably harder. Sections 3 and 4 
contain the requisite asymptotic analysis for the non-compact space . For small 
volumes, we use a Euclidean approximation (Section 3), and for large volumes we 
obtain the interesting result (Proposition l4.12"|) that 

lim F{iJW,w) = 27rln , 

which is positive if and only if > A = /A — 1 < .85. Along the line v = \w, 
F is decreasing in w, when w > 300. Moreover, once w > 150, dF/dv > 0. 
This asymptotic analysis reduces the problem in to a bounded domain. We 
are able to use a computer to show that the Hutchings function is positive on this 
bounded domain. In S^, which is a space of finite volume, our conjecture requires 
each region to contain at least 10% of the total volume of S^, and we are again 
examining a bounded domain. Next we use a new balancing argument (Proposition 
[SHU, which says that if F{^, ^) is positive, so is F{v,w) for v > 2w. This 
argument, together with symmetry (used only for S^) and a balancing argument 
due to Hutchings (Lemma [234}, reduces the size of the domain to check. We then 
break up the domain into small rectangles (and triangles). Finally, we decompose 
the Hutchings function as the difference of a concave function and an increasing 
function, so that it suffices to consider the values on (or beyond) the corners of 
the rectangles (see Lemma r5.ip . We find values of the parameters to locate such 
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points and verify that F{v,w) > 0, allowing a safe margin of error. If necessary, 
we subdivide the rectangles and repeat the process. 

Though using a computer introduces error, we bound this computation error 
and account for it. 

A more succinct version of this paper has been submitted to the Journal of 
Geometric Analysis [C 5] . 

1.2. Acknowledgments. The authors would like to thank Colin Adams for pos- 
ing this problem and for his insightful questions and helpful conversations. We 
would also like to thank David Puter, Matt Kudzin, and Pat McDonald for helpful 
conversations. Finally, we would like to thank Prank Morgan for his direction and 
advice. 

This research is the joint work of SMALL undergraduate research Geometry 
Group from '01, '02 and '03, and was completed by Neil Hoffman (Geometry Group 
'03) in his undergraduate thesis work [Ho]. We thank the National Science Founda- 
tion, Williams College, and the organizers of the SMALL REU for helping to make 
this research possible. 

2. The Hutchings function in and 

In and H^, given volume v, a round ball has least boundary area A{v) ([Sch], 
[M| p. 127]). Similarly, given volumes v,w, there is an area-minimizing double bub- 
ble that encloses and separates these two volumes; furthermore, this double bubble 
is comprised of smooth constant-mean-curvature hypersurfaces, except possibly for 
a set of measure zero ([Ml Thm 13.4, Remark before Prop 13.8], \CF\ Propostition 
2.3]), and it is symetric about some geodesic ([h1, Lemma 2.9, Remark 3.8]). Of 
course, in S^, we assume that v,w satisfy v + w < \S^\, where |S^| denotes the vol- 
ume of S^. Also, we put v,w,u to be w/|S'^|, w/|S^|, u/|S^|, respectively. Though 
the area of the area-minimizing double bubble is sometimes denoted A{v,w) we 
reserve this notation for the area of the standard double bubble. The following 
argument shows A{v, w) is well defined. 

Lemma 2.1 (Unique standard double bubble). Given volumes v and w: 

• Up to isometries of the space, there is a unique double bubble in 
consisting of three spherical caps meeting at 120° that encloses and 
separates a region of volume v and a region of volume w (whenever 
V + w < 1). 

• Up to isometries of the space, there is a unique double bubble in JrP 
consisting of two spheres (the outer caps) and a sphere, hyposphere, 
horosphere or geodesic plane (the inner cap) meeting at 120° that en- 
closes and separates a region of volume v and a region of volume w. 

Definition 2.2. The Hutchings function in is Fs^ : {{v, w) : v + vj < 1} ^ H, 
where 

V 

(1) Fs3{v, w) = As3{-) + As3{w) + As3{v + w) - As3{v,w). 

The definition of Fji^ is similar. 
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2.1. The Hutchings component bound. In this section, we recall for the reader 
the basic property of the Hutchings function, namely that it can be used to limit 
the number of components in an area-minimizing double bubble. 

We suppress the subscripts on F and A, since the results stated here apply to 
both S3 and H^. 



Proposition 2.3. If F{v, ui) is positive, then in an area-minimizing double 
bubble enclosing and separating volumes v and w, the region V is connected. 

Proof. This is a direct consequence of [CF\ Proposition 4.8] since the area of the 
standard double bubble has at least as much area as the area-minimizing double 
bubble. □ 

2.2. Area and volume formulas. This section recalls the formulas for area and 
volume of area-minimizing single bubbles (Remark 12. 4p . We have already shown 
the existence and uniqueness of a standard double bubble enclosing given volumes 
(Proposition l2.ip . followed by Propositions l2.61 12.7) 12.81 and 12. 91 which set up basic 
formulas needed to compute area and volume for standard double bubbles. 
Single bubbles. Here we present the standard results on volume and surface area of 
spheres in and H^. 

Remark 2.4. The formulas for surface area of a sphere of radius r in and H'^ 
are 

(2) As3 = 47rsin^r 

(3) Au^ = 47rsinh^r. 

The volume formulas for a ball of radius r in and are 

(4) Vs3 = 7r(2r-sin2r) 

(5) Vh3 = 7r(sinh2r - 2r). 

Lemma 2.5. The mean curvatures, ^ , of spheres of radius r in and 
are 2 cot r and 2 coth r respectively. 

Proof. Volume and area are related hy ^ = A. Thus, in S^, by Remark 12.41 



(6) 



dr 

dA dA/dr 



dV dV/dr 

(7) = A'{r)/A{r) 

, . Stt sin r cos r 

(8) = 

47r sm r 

(9) = 2cotr. 

The proof in H'^ is analogous. □ 

Since each standard double bubble in and is symmetric about a geodesic 
line, it is convenient to compute area and volume of a standard double bubble by 
considering the revolution of some generating curve consisting of three circular arcs 
meeting at 120°. Figure [3] represents an arbitrary generating curve labeled with 
relevant measurements. 



DOUBLE BUBBLES IN AND 



7 




Figure 3. Some important geometrical features of the generating 
curve for a standard double bubble. 

Proposition 2.6 (Area of Spherical Cap in S^). In the area of a spherical 
cap subtended by an angle 2(po is 

(10) / / sin'^ r sin0 (i^d^ = 27r sin'^ r(l — COS 0o)- 
Jo Jo 

Proposition 2.7 (Volume of Spherical Cap in S^). The volume of a spherical 
cap that IS subtended by an angle 2(po is 

(11) — 7r(tan^"'"(cos tanr-) — r cos 0o + (— 1 + cos ^o)(^ — cos r sinr)). 
Hyperbolic Space. 

Proposition 2.8 (Volume of Spherical Cap in H^). The volume of a spherical 
cap in a sphere of radius r subtended by an angle 2<po is 

(12) 7r(— r + tanh~'^(cos (po tanhr) — cos (po coshr sinhr + coshr sinhr))). 

Proposition 2.9 (Area of Spherical Cap in H'^). In the area of a spherical 
cap in a sphere of radius r subtended by an angle 2(po is 

(13) / / sinh^ r sin(/><i(/><i^ = 27r sinh^ r(l — COS 0o)- 
Jo Jo 

Now we can compute the volume and areas of standard double bubbles. 

For each bubble component, we compute the volume and area of the spherical 
cap which separates it from the exterior. To find the area of the double bubble, we 
simply add the area of the three spherical caps. 

To compute the volume of the bubble components we add the volume contained 
by the separating cap to that contained by the smaller of the two outer cap to 
obtain the volume of the smaller region, and subtract the volume contained by the 
separating cap from that contained by the larger of the two outer cap to obtain the 
volume of the larger region. (Of course, when the two outer caps are the same size, 
the separating cap contains no volume so there is nothing to add or subtract.) 

Examining Figure [3l we observe that the angle 0i has a different orientation 
(relative to and rs) depending on the relative size of the two bubbles. When 
01 is 90°, ^2 and are equal causing the volumes enclosed by both regions of the 
double bubble to be equal. Tracing back through the equations, we notice that this 
happens when 
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and 

(15) 1^ = 1/3 (inH3). 

We make adjustments to (pi accordingly, and call the adjusted angle 4>i- Using 
this angle in the formulas for area and volumes of spherical caps yields the desired 
result, i.e., if (pi > ^, the formulas return values corresponding to a cap smaller 
than a hemisphere, otherwise the values returned describe a cap that is larger than 
a hemisphere. 

2.3. Properties of Hutchings function. The concavity of A{v) is well-known 
in both and H^. Hutchings shows that the area of the area-minimizing double 
bubble enclosing volumes v and ti; is a concave function in both variables [HI Rmk 
3.8, Thm 3.9]. For our computations, the concavity of the standard double bubble 
is relevant. 

Proposition 2.10. A(v,w) in S" and H" n > 2 are strictly concave functions. 

Proof. By uniqueness of standard double bubble (Lem [27l|) . it is the only equilib- 
rium and hence the minimizer among spherical surfaces meeting along a sphere. 
Given Vq, uio, and the standard double bubble, consider the area Ai{v, w) of spheri- 
cal double bubbles for nearby values of v and w obtained by varying the curvatures 
of the two bubbles, but leaving the curvature of the interface unchanged and giving 
up the 120-degree-angle condition. Along any line in {v,ui) space, if v is strictly 
increasing, then the mean curvature Hi of the first spherical cap is strictly de- 
creasing, and similarly for w and H2. Hence A'l = Hidv/dt + H2dw/dt is strictly 
decreasing and Ai is strictly concave. Since A{v,w) < Ai{v,w), with equality at 
{vo,Wo), A{v,w) is strictly concave. □ 

Proposition 2.11. A(v,w) in is strictly increasing in v,w on the closed 

I I I I I I 

triangular domain defined by (0, 0), (■'-j^, ^g^), (0, -i-g^). Furthermore, A(v, w) 
is also increasing on the line segment with end points (0, ^^), (^^, ■^^)- 
Proof. Since A{v,w) is concave (Proposition I2.10|l and symmetric in both argu- 

I g3 I I g3 I 

ments, A{v,w) in S attains its maximum at (^,^). Since A{v, w) is also 

|g3| |g3| 

concave along lines through (3^, 3^) it is increasing on the triangular domain 

I I I S'^ I I S"^ I 

and the line segment with end points (0, ^j-*-), (3^, 3^)- ^ 

Proposition 2.12. For all v,w, Aii3{v,w) is strictly increasing in each vari- 
able. 

Proof. Aii3{v,ui) is strictly concave along any line, as shown in Lemma [2.101 in 
particular, along any line v = or w = Wq. Because it is positive for all v,ui > 0, 
it follows that it can never be non-increasing along any of these lines. □ 

2.4. Balancing. We call the following useful proposition "S-balancing" because it 
relies on the concavity of area of the standard double bubble, in contrast to Hutch- 
ings balancing (Lem [27l4|l which relies on the concavity of area of the minimizing 
double bubble. 

Proposition 2.13 (S-balancing). In R", H", or S", for n > 3, if w < v < 2w 
and F{^, ^) > 0, then F{v, w) > 0. 
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Proof. Consider an area pair v > ui and suppose that F{^^^y^, ^^^^) is positive, i.e., 
that 

(16) 2A {^-^j + A (^-^j + A{v + w)-2A > 0. 

Concavity and the symmetry of A implies that for any v not equal to w, 

'v + w V + 



> A{v,w) 



But then we have 



F{v, w) > 2A (-J + A{w) + A{v + w) - 2A ( , 
So it suffices to show that 

'v + w v + w 



2A (^^ j + A{w) + A{v + w)-2A 
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, / V + w\ , I V + w\ ,, . , f V + W V + w 

> 2A[^-^yA[^^yAiv + w)-2A[^-^,^ 

i.e. that 

, , . . f V + w\ , ( V + w 

(17) 2Ai^-)+A{w)>2A[-^yA[-^ 

Consider three spheres, two enclosing volume a;/2 < y and one enclosing volume 
y. Initially, increasing x and decreasing y while keeping x + y constant will increase 
the total surface area (we know this because 2A{xl2) + A{y{xy) is a concave func- 
tion). Indeed, area will increase until a;/2 = y, that is, until all three spheres have 
equal volume and any further transition would be wnbalancing. In particular, if 
w < V < 2w and we put x = y = {v + iy)/2, then increasing x until it equals v 
while decreasing y until it equals w will increase the total surface area of our three 
spheres. This shows (Equation [TT]). 

Hence F{v, w) > for w < v < 2w. 

□ 

Lemma 2.14 (Hutchings Balancing). In R", H", or S", for n > 3, if v > 2w 

( or in if v > 2u ), then the region of volume v is connected. 

Proof. See [H, Remark 3.8, Theorem 3.5, Corollary 3.9]. □ 

2.5. Permutation. Permutation utilizes symmetries that exist with the Hutchings 
function when applied in compact spaces. Though this argument is applied in S^, 
it can be used in any space with finite volume that has a Hutchings function. 

Proposition 2.15 (Permutation). F{v,w) = F{v,u) where u is the volume of 
the complement. 

Proof. First notice that A{v, w) = A{v, u) as it is only a matter of labeling U or W 
the exterior. Now, A{v + w) = A{u) and A{v + u) = A{w), since U = {VUW)^, and 
similarly, W = {VUU)^. Thus, F{v, w) = 2A{^) + A{w) + A{v + w) - 2A{v, w) = 
2A{^) + A{w) + A{u) - 2A{v, w) = F{v, u), as desired. □ 
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Remark 2.16. S-balancing (Prop [27T3|l . Hutchings balancing (Lem l2.14|l , and per- 
mutation (Prop irTS]) yield a complete answer to a question of Cotton and Freeman 
[CF\ Open Question 1]. In particular, the larger region of an area- minimizing dou- 
ble bubble in is always connected. If v > 2w, this follows from Hutchings 
balancing (Lem [27l4|l . Otherwise, connectedness follows from the positivity of the 
Hutchings function F{v, w) (Prop [273|) . This in turn follows from S-balancing (Prop 
I2.13P and the positivity of the Hutchings function F{v,v) for equal volumes ^CFl 
Prop 5.11, 5.14,5.19]. In S^, Hutchings balancing (Lem [27l4|) shows connectedness 
when V > 2vj oi v > 2u. By Proposition l2.31 it suffices to show that F{v, w) > if 
w <v <2w oiu <v <2u. If w < .4, F{v, v) > Q ([CP, Prop 5.5, 5.8]). S-balancing 
(Prop [2713]) shows that F{v,w) > ii w < v < 2w. Finally, permutation (Prop 
[27T51I shows F{v, w) > ii u < V <2u. 

3. Positivity of the Hutchings function for small volumes in 
3.1. Properties of the Hutchings function in H^. 

Theorem 3.1 (Double Bubble Theorem in [HMRR] ). The unique least-area 
enclosure of prescribed volumes v,ui in R^ is a standard double bubble. 

Remark 3.2. For small prescribed volumes in H'^, the space around any standard 
double bubble is nearly flat and looks like a portion of R^. By estimating the 
amount of distortion carefully, we can use information about the Hutchings function 
in R^ to prove that the Hutchings function in is positive for certain prescribed 
volumes near 0. In the following proof, A is a constant which captures the amount 
of distortion, and r represents the radius of a small ball in which is only slightly 
affected by the curvature of . 

Lemma 3.3. For any pair of prescribed volumes less than (^ + ^^)(^:f^^)^) 
some minimizing double bubble in enclosing those volumes fits inside a 
ball of radius r. 

Proof. Without loss of generality assume that v < w and flx w. The exterior 
spherical cap enclosing the region of volume w has a radius, r^. Consider standard 
double bubbles where the radius of this spherical cap is flxed at rw In this case, 
the equal volumes double bubble has the greatest diameter. However, the region 
enclosed by the spherical cap with radius has volume less than w. By increasing 
the radius until this region has volume w, the diameter increases. Thus for any 
non-equal volume double bubble enclosing volumes v and ui, where v < w, there 
will be an equal volume double bubble enclosing v' and w with a bigger diameter. 

If is the radius of the spherical cap that forms the exterior of W, then the 
diameter of the equal volume double bubble is 2rw + ryj\/3. So, this double bubble 
will flt inside a ball of radius r > r^ ^■^'^^^ . 

We also know that in R'^ a spherical cap of radius ro subtended by an angle 2<p 
has volume: 



(18) 



V{r) = -^i'^o)^ + Trro^cos<p. 
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Sincere [0,f], 



3 27r ^ 7r-\/3^ 



(19) w<{r^y{—+ 2 
Substituting "^^ r for in the above equation gives the desired result. 

2+ V 3 

□ 

Lemma 3.4. If .84w < v < w, then 

2Ar^{^) + Ar^{w) + Ars{v + w) > (2.02676)^R3(i;,w) 

holds. 

Proof. (Modification of [HMRR| Proposition 6.2]) By rescaling, we may assume 
that V = 1 — w and that w £ [|, j^g]. Because Aji3{v,w) is concave [HI Theorem 
3.2], we have 

^R3(l -W,W)< ^R3(i, 1) = (2-t • 3)^R3(1). 

Hence, it suffices to prove that 

2AR3{i^) + An3{w) + Ar3 > 2.02676(2-1 •3)^r3(1), 

or equivalently (dividing by j4r3(1) on both sides) that 

23(1 -iu)i +wi +1> 2.02676(2-t • 3). 

At w = I and w = Ykg, the left hand side of the above inequality is bigger than 
2.4236 and 2.412965, respectively. The right hand side of the above inequality is 
less than 2.412966. Hence, the inequality holds for ui = ^ and w = and 
because the left hand side is concave in w, the inequality holds for all ui £ [|, 
as desired. 

□ 

Proposition 3.5. IfO<v,w< .002743 and .84w <v <w, then 

(20) 2^H3(^) + -AH^iw) + j4h3(w + w) - 2Aii3{v,w) > 0. 



Proof. By Lemma [3731 we can fit a double bubble in R'^ with prescribed volumes 
v,w < .002743 in a ball B of radius r < .1547. 

For ease of notation let A = ^i^^ < 1.003994. 

r 

Notice that if this ball can enclose a standard double bubble of volumes {X~^v, X~^w) 
then it can enclose spheres of volumes: 

{A"^^, X~^w, X~^{v + w)}. We may estimate the surface area of any sphere in 
which fits inside B. Suppose that a sphere enclosing volume X~^a fits inside B; 
in particular, this holds for a £ {^,ui,{v + w)}. The single bubble in encloses 
a volume a' < a and has surface area at least j4R3(A-'^a). Thus, 

(21) ^H3(a) > ^R3(A-2a). 
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Therefore, 

2^H3 (|) + ^H3 (w) + Ah^ {v + w) 

> 2Ar3{X-^v/2) + Ar3{X-^w) + Ar3{X-^{v + w)). 

By Lemma [331 

2Ar3{^-) + Ans{X-^w) + Ans{X-^{v + w)) 

> (2.02676) Ar3{X-^v,X-^w). 
Scaling of double bubbles in tells us 

(22) Ars{X-\,X-^w) = (X)-^/^Ars{v,w). 

Because B contains some minimizing double bubble Si in enclosing volumes 
v,w, we can estimate the surface area of the image of Si in H^. This double bubble 
S2 encloses volumes V2 > v and W2 > tw, with surface area less than or equal to 
A-Ri{'U2,W2)X^ . Therefore, 

(23) ^R3(w,w) > A-2^H3 (1^2,^2). 

By Proposition [57121 -^ti^i^^ v) is increasing in each variable, implying that 

(24) ^H3 (1^2,^2) > Aiiz{v,w). 

Combining inequalities ([22]), ([22]), ([23|, and ([24} with equality ([22]) gives 

(25) 2Aii3{^) + An^iw) + An^lv + w) > {2.0676){X)-^°/^Ah3{v,w). 
Note that (2.02676)A-i°/3 > 2, so 

(26) 2^H3(^) + Aii3{w) + Aii3{v + w) > 2Aii3{v,w) 
as desired. 

□ 

Remark 3.6 (Notes on S^). The above argument works for as well, although 
we do not need it because we are just concerned with volumes greater than 10% of 
the total volume of S'^. The proof uses a different distortion factor, A = ^^j^. 
Otherwise the proof follows with very little adjustment. Using Lemma 13. 4| the 
Hutchings Function is positive for small volumes in namely v,w < 0.002738 and 
.84w < V < w. 

4. POSITIVITY OF THE HUTCHINGS FUNCTION FOR LARGE VOLUMES IN 

4.1. Introduction. Theorem 14.231 shows that for sufficiently large v and w, the 
Hutchings function in is positive whenever tp = v/w > X and w > 300. 
Throughout this section we set A = 1. The proof of Theorem 14.231 uses 

three ancillary results: 

Proposition [4712] . For a fixed ratio v/w, the limit of the Hutchings function as 
w ^ 00 is nonnegative if and only if ip > X < .85. 

Proposition [4720] • The Hutchings function strictly decreases if one travels out- 
ward along the line v = Xw as long as w > 300. 

Lemma [4.221 . The partial derivative with respect to v of the Hutchings function 
is positive if w > 150, v > Xw. 
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Together, these results show that the Hutchings function is positive for tp yX and 
w > 300. In order to prove Proposition l4.12| we must describe the Hmits of various 
quantities for standard double bubbles enclosing volumes v^w as v,w grow large. 
Proving Proposition l4.2Ql and Lemma [4. 221 require a closer examination of how fast 
these quantities converge to their limits. So, the results about the limits appear 
alongside various inexact numerical estimates. 

4.2. Preliminary material: spheres in H^. If r = r{v) is the radius of a sphere 
with volume v, then 

(27) 2w/7r = e^^ -4r < e^'' 
so that 

1 , 2v 

(28) '">:^ln-. 

Lemma 4.1. If r = r{v) is the radius of a sphere with volume v, then 

(29) lim (r - - In — ] = 0. 

v—Kx> y 2 TT y 

Proof. By manipulation of the volume formula, 

(30) ^(e2(r-iln2.)_^^^g_2. 

TT 

Notice for any e > we can find a sufficiently large r, such that, 

(31) e"2r < 2(sinh2r - 2r){e^' - 1). 

Hence, 2Ji(e2(r-i in 2j^) _ jg iggg tji^n f{e^' - 1) and 2(r - |ln^) < 2e, as 
needed. □ 

Remark 4.2. When v > 15QA > -7r(sinh(4.4) -4.4), we have r > 2.2 and Inequality 
([3T|I holds for e = 0.06. Thus, r < | In ^ + .06 for such v. We obtain the following 
numerical estimate: if > 150A < 127.09, then r < Mn — + .06. 

Lemma 4.3. As v approaches infinity, A{v) tends to 2i; + 27rlnt; — 27r(l — 
ln(7r/2)). 

Proof. Let r = r{v) denote the radius of a sphere with volume v. Since 

(32) r = -In ( — + 6^2'' +4r 

2 \ TT 

we see that, 

(33) A{v) = 2v-2-K + 27re~^'' + 47rr. 



By Lemma mil 

\ 2v . 1 1 TT 

(34) lim r In — = lim r InwH — In — = 0. 

^ ' v^co 2 TT v^oo 2 2 2 

Also, limt;^oo e^^'" = 0. The desired result follows easily. □ 

Lemma 4.4. The curvature of a sphere with volume v and radius r is 

, , 27r 

(35) 2 + 



t; + 27rr + f e-^r- _ 7r/2' 
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Proof. The curvature equals 

-I- e"'' 4e~'' 4 

(36) 2cothr = 2 = 2 + = 2 + -^ . 

Also, 2v/Tr = e^'' — e~^^ — At. The desired result follows easily. □ 
Lemma 4.5. A!{x) < 2 + -^^2^ — 5. 

Proof. Prom ([28]) we have the following bound on the radius r of a sphere containing 
volume x: 

(37) r > ^ln(2a;/7r) 

(38) = ^lna; + (l/2)ln(2/7r) 

(39) > ^Ina; - 3/27r+ 1/4. 
By Lemma Sm 

-^'(^) = ^+ . + 2.r + (J/V---./2 - 

Also, e~^'^ > 0. Applying these inequalities in the above expression for A'{x) yields 
the desired result. □ 

Lemma 4.6. If x > 150A then A'{x) > 2 + ^+^i^^J_i.04i • 
Proof. By Lemma SHI 



(41) A'{v) =2+ 



v + 27rr + f e-2'- - f 



The denominator is a sum of three terms: 

-2r 

' 2 2 

Because e~^^ < 1, the third term is less than 0. We now approximate the second 
term. By Remark |4.2| we have r < |ln^ + 0.06 because v > 150A. Hence, the 
second term satisfies 

2v 2 

(42) 2-7rr < Trln h 0.127r = Trlnw + Trln- + 0.127r < -/rlnw - 1.041. 

TT TT 

With these approximations, we find 

.1, N 2'K 

(43) A\v) > 2 + 



w + Trlnw - 1.041 

□ 

Corollary 4.7. If x > 300A t/ien yl'(a; + 3) > 2 + ^^^p^^- 
Proof By Lemma ^'(i) > 2 + ,+^i,';_i.04i , if a: > 150A. 

Because e^^'^ < 1, the third term is less than 0. We now approximate the second 

term. By Remark |4.2| we have r < |ln^ + 0.06 because v > 150A. Hence, the 

second term satisfies 

2v 2 

(44) 2-7rr < Trln h 0.127r = Trlnw + Trln- + 0.12Tr < Trlnv - 1.041. 

TT Tr 
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With these approximations, we find 

2-K 



(45) A'iv) > 2 + 
^ ' ^ ' V + -Klnv - 1.041 

For X > 300 A, we have 

(46) TT In ^-i-^ - 1.04K -1.004. 

X 

Hence, 

, , ^ 27r 27r 

(47) A'(x + 3)>2 + ; >2 + ; . 

^ ^ ^ ^ a; + 3 + Trlna; - 1.004 a; + 7rlna; + 2 

□ 

4.3. The Hutchings function has a positive Hmit. 

Lemma 4.8. Given a standard double bubble, the radius of the circular inter- 
face where the three caps meet is less than cosh^"'(2). 

Proof. Consider the disc whose boundary is the circular interface. Let y be the 
radius of this disk. It meets the bigger of the outer caps of the double bubble at 
an angle a > 27r/3. The radii of this outer cap and the disk come together at an 
angle of /3 > 7r/6 along the circular interface. Let 2<po be the angle that subtends 
this cap. By the formula for the area of a spherical cap in (Proposition 12. 9p 
and the hyperbolic laws of cosines and sines [T] and [St] , 

(48) cos 00 = — cosh J/ sin /3 

• , sinhj/ 

(49) sinhr' = ^ — 

sm cpo 

the surface area of the outer cap is 

, . 27r sinh^ y 

(50) 



1 — cosh y sin yS 

The surface are of this cap is positive and finite, implying that coshy < csc/3 < 
csc(7r/6) = 2. Because t cosht is increasing for t > 0, y < cosh^"'(2). 

□ 

Lemma 4.9. Given a spherical cap subtended by an angle of2(p on a sphere of 
radius r, let 6 be the angle the sphere makes with the disk of radius y bounding 
the cap from below. The volume of this cap can be expressed as, 

. . ( sinhwsin^ , -i f sinhj/sin^ 

(51) TT ; - tanh 

\sechy -\- cos 6 \coshj/ + cos6 

Proof. This follows from using the identities: 

tanh(a + 6) = ^ifj+^tl^^, where (a + 6) = (-r + tanh" ^ (tanh r- cos 0)), 
1 — tanh'^ X = sech'^o; and cosh 6 = ^15^ 

cos 7 ' 

cosh c = cosh y cosh 6, 
cos 7 = sin cosh fe, 
tanh r cos 7 = tanh y, cosh y sin 7 = cos <p, 



16 



J. CORNELI ET AL. 



sinh y = sin cp sinh r and cos 7 = i^HkM ^ and 



sm — 



tanh r ' 

sinh r ' 

□ 



Lemma 4.10. Given a standard double bubble enclosing volumes v, w, suppose 
that 6 is the angle between the separating cap and the disc with the same 
boundary, and let y be the radius of the disc. As v,w grow large, 6 approaches 
and y approaches cosh" ''(2). 

Proof. By Lemma [4.81 y < cosh^"'(2). The two caps of the region containing the 
disk meet the given disc of radius y at angles 6 £ [0, 7r/3) and a = 2Tr/3 — 9. By 
volume formulas, the volume x of this region is the sum of 

, . / sinh w sin 6 , _■, ( sinh j/ sin ^ 

(52) — — — — - tanh 



\ sech y + cos 6 \ cosh y + cos B 

and 

/ sinh J/ sin Q / sinh y sin a 

(53) TT — ; tanh ; 

\sechy + cosa \coshj/ + cosa 

Notice that coshy + cos^ > 1 H — 1 = 0, and similarly coshy + cos a > 0. Hence 
the "tanh"^" terms are positive, implying that 

, . / sinh y sin a sinh y sin ^ 

(54) a: < TT -— — —— + 



\sechj/ + cosa sechj/ + cos6^ 

Because sechj/ > 1/2 and cos6,cosa > —1/2, the denominators of the above 
expression are positive. Hence, 

'sinh(cosh"^(2)) sin(7r/2) sinh(cosh"^(2)) sin(7r/2) 
sechj/ + cos0 sech y + cos a 



< 



y sech y + cos Q sech y + cos a 
27rV3 



sechy + cos(27r/3 — 6) 

As V, w approach infinity, x > m.m{v, w} must approach infinity as well, implying 
that the positive quantity sechy + cos(27r/3 — 6) must approach - that is, y must 
approach cosh" (2), and cos(27r/3 — 9) must approach — |. Because 27r/3 — ^ £ 
(7r/3, 27r/3], 6 must approach 0. In particular, when v,w > Xq = 300A, we have 

^ 27rN/3 27rx/3 



sech y + cos(27r/3 - i9) 1/2 + cos(27r/3 - fi^) 
or 

Stt /3 

(56) cos(27r/3 - < — 1/2 < -.457. 

Xq 

For e e [0, 7r/3), this is true only if 6^ < 1/20. □ 

We have the following numerical estimate: ii v,w > 300A ss 254.18, then 6 < 
1/20. 
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Lemma 4.11. As v,w grow large, the surface area of a standard double bubble 
enclosing volumes v,w approaches 

(57) Aus{v + Voo) + Aii3{w + Voo) - '2aoo + Coo 
where v^o = 7r(3/2 — ln2), aoo = Stt, and Coo = 27r. 

Proof. Given a standard double bubble enclosing volumes v,w, consider the region 
enclosing volume v. Its outer cap lies on a sphere; let the volume of this sphere 
be w + vi]w and let the surface area of the sphere exceed that of the outer cap by 
ai]lu. Define vi% and a[^i similarly with respect to the region enclosing volume 
w. Finally, let Cy^w denote the surface area of the separating cap. Then the surface 
area of the standard double bubble is 

(58) {v + vi]i) + Aus {w + vi%) - a(i) - -(2) 



v.w 



Suppose that 6 is the angle between the separating cap and the disc with the same 
boundary, and let y be the radius of the disc. By Lemma [4.101 as v and w grow 
large, 6 — > and J/ — > j/o = cosh^"'(2). Using formulas for volume and surface area, 
one can easily check that 

lim v^^^ = lim w^^^ 

/ sinh j/o sin(7r/3) t -i f sinh j/o sin(7r/3) 

TT — 7 — — — tanh y — — 

\sechy + cos(7r/3) \cos y + cos(7r/3) 



lim all' = lim aP' 



27r sinh^ j/o 



and 



lim Cy ■ 



1 + cos(7r / 3) cosh j/o 



27r sinh^ j/o 



Hence, 



v,w^oo ' 1 + COS cosh 2/0 

— Coo . 



\A{v + v'^^l) - A{v + Voo)\ < A'{v)\vi'X 



(2) 

approaches as v,w grow large. So too does \A{w + Vy^i,) — A{w + Vao)\- From 
these computed limits, it follows that the surface area of the standard double bubble 
enclosing volumes v,w - as calculated in Equation [58] - approaches A{v + Wqo) + 
A{w + Voo) — 2aoo + Coo as v, w get large. □ 

Proposition 4.12. For each fixed ip > 0, 

4(t/' + 1) 



lim F{tpw,w) = 27rln- 



=2 



(Note that this limit is nonnegative if and only if Tp ^ ^ which is less than 
0.85.; 
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Proof. We write v = v{w) = \w. Let Wqo = 7r(3/2 — lii2), = Stt, and Cqo = 27r. 
By Lemma [4.111 the given limit equals 

(59) lim 2A{^) + A{w) + A{v + w) - 2{A{v + v^) + A{w + Voo) - 2a^ + Coo) 



provided that the latter limit exists. By Lemma [4731 

(60) lim A{x) = 2x + 27rlna; - 27r(l - ln(7r/2)) 

a:— >oo 

fori e {^,w,v + w,2v + Voo,vj + Voo}- Substituting these expressions into [59l and 
simplifying gives 

lim (27r(21n(-) + Inw + ln{v + w) - 2ln{v + Voo) - 2ln{w + Voo)) 

w^oo 2 

-8voo + 4aoo - 2coo) 
provided that this new limit exists. Indeed, this new limit does exist and equals 

hm 27rln -— \^ ] - Sv^o + 4aoo - 2coo 

f V + w\ 

= lim (27rln — - 8voo + 4aoo - 2coo) 

= 27rln ( ^^^^ ) - 8 • 7r(3/2 - ln2) +4 • Stt - 2 • 27r 



2-7rln 



4 



Hence the limit in ([59]) evaluates to 27r ln(^^^^^^), which is only positive when 
ip ^ ^, finishing the proof. □ 

4.4. The Hutchings function is decreasing along a line for large volumes. 

Lemma 4.13. Given a standard double bubble enclosing volumes v,w, the 
outer cap of the region of volume v lies on a sphere, say with volume v + Vi 
(reap, the outer cap of the region of volume w lies on a sphere with volume 
w + Wi). If v,w > 300X then Vi,wi < 3. 

Proof. Let vol(j/, a) denote the volume of a spherical cap with associated disc of 
radius y, where the cap and disc meet at an angle [0,7r]. This volume is increasing 
in both y and a. Without loss of generality, assume that v < w. The three caps of 
the standard double bubble meet along a circular interface with radius j/o; let the 
disc with this circle as its boundary meet the separating cap of the double bubble 
at angle 6. By Lemma yo < cosh"^(2). By Lemma OOl 6 < 1/20. Both 
Vi = vol(2/o , tt/S + 6) — vol(yo , ^) and 11)1= vol(j/o , 7r/3 — 6) + vol(j/o , ^) are less than 
vol(2/o , 7r/3 + 6) + vol(2/o ,6), which in turn is less than vol(cosh~ (2) , 7r/3 + 1 /20) + 
vol(cosh"^(2),l/20) < 3. □ 

We now prove some algebraic lemmas that will be used in the proof of Proposition 
14.201 and later in Proposition l4.21l 

Lemma 4.14. If ttItlx + a > 0, then ^ ^ > ^ - 

' ^ ) x-\-7r \n.x-\-a x x^ 
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Proof. The left hand side equals 

, . 1 TT In a; + a 

X a;(a; + Trlna; + a) ' 

which is greater than ^ - 2Jn|+a 

Corollary 4.15. For positive numbers jj,, w with In^u -I- a > 
, . u 1 l/7rlnw Trinu a 

(62) , . > - - — + 



juw + TT In(juw) + a w w-^ \ jj, jj, jj, 

Proof. Simply apply Lemma [4. 141 with x = fiui. 



Lemma 4.16. If x > 150A, then 



1 1 TT In a; 

< - 



i + Trlna; — 3 x l.lx'^ 

Proof. For such x, we have l.li >a; + 7rlna; — 3>0. Hence, 

1 1 Trim — 3 iTrlna; — 3 

a; + 7rlna; — 3 x a;(a; + 7rlna; — 3) x l.lx'^ 

as desired. 

Corollary 4.17. For positive jj,, w , with iiw > 150 

A* < 1 _ J_ TTlnw _^ TTlnyU 3_ 

tXW + 'K\D.{fXw) — 2i W W^ \ l.ljj, 1.1/i 1.1/i 

Proof. Apply Lemma [4. 161 with x = fiw. 

Lemma 4.18. For all w > 300, 

2^^ + ^--+.ln..-2j< 
27rlnw 27rln(|) 6 Trlnw 3 

1 lAl. 1 L 

l.lf l.lf l.lf 1.1 1.1 

Trlnw ^Trln(A + l) 3 



1.1(A + 1) 1.1(A + 1) 1.1(A + 1)' 
Proof For all w > 300, 



^(2^-2-^''^/^ ln(A + l) 



A l.lA/2 l.l(A + i; 

+ ''Til7^ + rr + rulTT)-A-^i< 

2 1 1 2 , , 

H 1 TT r - - - 2 Trlnw. 



l.lA/2 1.1 1.1(A+1) A 
Simple algebraic manipulation gives the desired result 
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Lemma 4.19. If w > 300, then 



2A 

4 + 4A + 27r — - — ^ h 



Xw + TT hi{\w) + 2 w + TT In(iy) + 2 
> 4 + 4A + 27r ( -— + ^ 



+ - 



Af + 7rln(Af)-3 w + '!rln{w)-3 
A + 1 



(A+l)u; + 7rln((A + l)iy)-3 
Proof. By Corollary SUll 



2A 

+ 



Xw + TT In(Aw) +2 w + TT ln(ii/) + 2 

f2 1 /Trlnw 7rln(A) 2A 1 , , 

\w w' \ X X X J vj^ 

4 1 / 27rlnu; 27rln(f) 6 Trlnu; 3 
up y l.lA 1.1 lA 

^ Trlnw ^7rln(A + l) 



1.1(A+1) 1.1(A + 1) l.l(A + i; 
(by Lemma OSjl 

A 1 (1 + A) 

> 1 h 



2|i + 7rinM_3 ' u;4-7rlnw-3 ' (A + l)u; + tt ln((A + l)ty) - 3 

(by Corollary I4.17p . Multiplying both sides of the inequality by 27r and adding 
4(1 + A) yields the desired result. □ 

Proposition 4.20. If w > 300, then 

-^F{Xw,w) < 0. 
aw 

Proof. Consider a standard double bubble enclosing volumes Aw and w. Let Vi,wi, 
be the volumes needed to complete the two outer caps of the double bubble which 
enclose volumes Xw and w. Then, 

dA 

(65) (^^i ^) = A^'(Ai(; + vi) + A'{w + wi). 
dw 

The bound on A'{x) given in Lemma [4751 tells us that 

(66) (2 + p. r -]>A'{^iw). 

^ ' V /iw + TT ln(/iiy) - 3 7 ^ ' 

Hence, 

(68) 2 + 27r ( ) > ^'(w). 
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(69) 2 + 27r(- r ^^-^ ] > A'((X + l)w). 

By Corollary SlTl 

2-K 



(70) A^'(Au; + 3)+^'(u; + 3)>A 2 



Xw + '7rln(Au;) + 2 



27r 

+ 2 + 



w + TT In w + 2 



Since Xw and w are greater than 300A, we can apply Lemma [4.131 to see that 
'^1,^1 < 3. Hence, 

(71) A^'(Au; + vi) + ^{w + Wi) > XA'{Xw + 3) + A'{w + 3). 
Therefore, 

(72) XA' {Xw + vi) + A' {w + wi) > 

27r \ / 27r 

A 2 + — + 2 + 



Xw + TT In Aw + 2/ \ ti; + 7rlnw + 2 

Combining inequalities [72l \67\ HH [69] and Proposition I4.19| we can see for all 

Xw,w > 300A, 

(73) 2XA'{Xw + vi) + 2A'{w + wi) > XA' ( — ] + A'{w) + A'{Xw + w). 



2 



Thus, 



(74) 2^A{Xw, w) > A^' f — ) + A'{w) + A'{Xw + w) 

aw \ 2 J 

as desired. □ 

4.5. The Hutchings function is increasing in v for large volumes. 
Proposition 4.21. For v > Xw, w > 150, 

^ I }^ ^ 2 ^ 

^ ' f +7r In f - 1.041 w + w + 7rln(v + !(;)- 1.041 vkIhv - 3 

Proof. Since v > Xw, w > 150, 



(76) y > 4(w + w)(7r In ^ - 1.041) 
and 

V -\- w 

(77) — - — >7rln(w + iy). 
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Hence, 

, . . , TT In ^ - 1.041 „ 

(78) v^{v + w)> {v + wY 2_ h w TT ln(v + lu) 

4 

, , . , TT In V — 3 
- 1.041^2 -{v + wf 



1.1 

So, 

2 7rlnf-1.041 1 tt ln(t; + to) - 1.041 1 7rlnw-3 
(I? "^^rr^ {v + wf ^'^^v l.lt;2 J 

Applying Lemmas 14. 141 and 14. 161 to the above equation yields, 
, , 1 1 1 

f + 7rln| - 1.041 w + w + 7rln(w + iti) - 1.041 w + Trlnw-S 

as desired. 



Lemma 4.22. For any fixed ui > 150 and v > \w , 

, , dF 

Proof. Differentiating the Hutchings function with respect to v yields 

(82) ^ = A'C^) + A'{v + w)- 2^A(v, w) > A'(^) + A'(v + w) - 2A'{v) 



□ 



By Lemma [4!6l and Lemma |475| 

, , dF 1 
(S3) 



dv I + 7rln^ - 1.041 

1 , 1 



+ 



V + w + TT hi{v + w) — 1.041 V + TT In V — 3 



By Proposition l4.21| 



, , dF 
(84) — > 0, 



as desired. 



4.6. Conclusion. 



dv 

□ 



Theorem 4.23 (Hutchings function is positive for large volumes). For all v < w 
such that vj > 300 and v > .85w ( or indeed v > Xw .841w ), the Hutchings 
function F{v, w) = 2A{^) + A{ui) + A{v + w) — 2A{v, w) is positive. 

Proof. Consider {v,w) = {vo,Wo), where A = Vq/wq > .85 > A and Wq > 300. By 
Proposition 14.221 the derivative of F with respect to its first argument is positive 
for the volume pairs under consideration, so F(wo,Wo) > F{Xwo,Wo). 

By Proposition l4.20l F(Xw. w) decreases as w increases and by Proposition l4.121 

i^(Au;o,u;o) > lim F{Xw,w) = 0. 
Hence, F{vo,vJo) > 0, as claimed. □ 
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Remark 4.24. The ratio A = — 1 is sharp in the following sense: for tp < X, 
there exist arbitrarily large v, w with v/w = ip such that the Hutchings function is 
negative. This is a result of Proposition l4.12l 

5. The computer proof of the positivity of the Hutchings function in 

S3 and H3. 

5.1. Introduction. Rigorously showing positivity of the Hutchings function is 
equivalent to showing that g{v, w) = 2A{^) + A{w) + A{v + w) is strictly greater 
than h{v,w) = 2A{v,w). The functions g and h have very similar properties in 
and H^, but the proofs vary slightly because the area of single bubbles and 
standard double bubbles are not increasing for all volumes in and the code 
has to efficiently calculate the curvatures of double bubbles. 

The computer proof relies on the facts that g is concave and h is increasing 
on relevant domains (Propositions 12.111 and I2.12P in both spaces and that g is 
increasing on all of H^. If g and h were each concave increasing functions of one 
variable, we could show that g is bigger than h on the interval [t;, w + 6], by showing 
g{v) > h{v + a) and g{v + a) > h{v + 6) (See Figure[4|. 

The computational lemma (Lemma IS.ip extends this idea into two dimensions 
showing that 5 > ft on a polygonal domain in the vw-plane if the minimum value 
of g on the vertex set is bigger than the maximum value of h. 

Using a computer does introduce error, which must be accounted for. A com- 
puter uses approximations for g and h, which we will call gcomp and hcomp- How- 
ever if g is underestimated and h is overestimated then the computer test will show 
gcomp > hcomp Only if ^ > ft. 

The following lemma shows that checking a finite number of points can give a 
global inequality. 

Lemma 5.1. Given a polygonal domain D with vertices identified by ordered 
pairs {vi,Wi), a concave function g{v,w), and another function h{v,ui) that 
obtains its maximum at {vo,Wo), if min{g{vi,Wi)) > h{vo,Wo) then g > h on 
D. 

Proof. Since g is concave it has a minimum on Z? at a vertex of the polygon. Thus, 
min{g{vi, Wi)) is a lower bound for g on D and h{vQ, Wq) is an upper bound for ft 
on D. Hence, g > h follows directly. □ 

The computer proof in S'^. In S^, rectangular regions and triangular regions are 
considered. There is additional computational complexity in because A{v) and 
A{v,w) are not increasing for all volume pairs {v,w) in S^. For rectangles where 
h{v,w) is strictly increasing in a neighborhood of the rectangle, the area of the 
standard double bubble is over-approximated by using radii of the outer spherical 
caps corresponding to volumes V4 and W4 bigger than v and w respectfully (See 
Figure p. 

For the area of the single bubble each part of g{v,w) is under-approximated. 
j4(|) and A{w) are always under-approximated by taking radii that correspond to 
spheres enclosing volumes less than | and w. Evaluating gi at (t;i,wi) in Fig- 
ures [5] and [7| gives an lower bound on gi{v, w) on the whole polygon. A{v + w) is 
under-approximated by considering the minimum of the value of j4 at a point to 
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the lower left and a point to the upper right of the rectangle. This is the minimum 
of p2(^'2,W2) and 52(^3, tfa) in Figures [5] and [3 

If g{v,w) < h{v -\- a2,'w + 62) (where w + 02 and + 62 correspond to Figure 
[5]), then the inequality is checked on the four rectangles with lower left corners of 
{v,w), {v + ai,w), {v,w + 61) and {v + ai,w + 61). 

Next, we consider the entire region where numerical analysis is necessary (Figure 
[6]). The triangular domain is first broken into parts: the right triangle {v,ui) = 
(.1,.1),(.1,|), (|,|) and the triangle {v,w) = (.1,1), (1,1) and (.1,.45). The lat- 
ter, upper triangle is difficult to deal with because h is increasing up to the hy- 
potenuse and decreasing above the hypotenuse. We can assure that we land exactly 
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on this line by setting the two parameters ri and to be equal. Getting an over- 
approximation for area on this line just involves overestimating ri. In addition, 
j4(|,|) must be hard coded into the program. If we want to over-approximate 
h{v, w) we move the point {v, ui) southeast along the line where w = u to (^4, UI4) 
and use h{v4, W4) as an over-approximation (See Figure[7|. For g, we still just take 
the appropriate corners of the rectangle containing the triangle. If not, then we 
break the triangle into a rectangle and two smaller right triangles by connecting a 
point on the hypotenuse to the two legs of the triangle via perpendicular lines. The 
desired inequality is then shown to be true on the rectangle and the two smaller 
triangles. Similarly for the rectangle if the desired inequality can not be shown 
on the large rectangle, the rectangle is split into four smaller rectangles and the 
inequality is shown to hold true on them (See Figure [5]). Both subdivisions repeat 
until the desired inequality can be shown on each subdivided rectangle or triangle. 
In the case where the upper right corner a rectangle is on the line where ui = u, 
we over-approximate v,w using the same methods as the triangle. 

The lower triangle can be replaced by the rectangle containing it, which can 
then be subdivided into smaller rectangles by the argument above. 

Mathematica uses computer algebra (with infinite precision variables). When 
checking inequalities, it reports true, false, or null if it cannot tell for sure |Wj. 

The computer proof in H^. In H'^, rectangular regions that tile the area between 
two lines through the origin are considered (See Figures [8l and [TO|l . The function g 
is strictly increasing, so g{vo,uio) is less than any value of g{v,ui) on a rectangular 
region with Vq < v and Wq < w. We can assure that these two conditions hold if 
we take a guess at the curvature that we know to be too low (corresponding to a 
volume that is too high) . Then, we increase the curvature (decreasing the volume) 
of the sphere by multiplying it by some constant. This process is done first with a 
relatively large constant until the curvature corresponds to volume that is too low; 
then the curvature is divided by the constant and the curvature is increased by a 
smaller constant until it corresponds to a a volume that is less then the volume 
being approximated. Using this curvature for our surface area calculations under- 
approximates the surface area of single bubbles. 

For A{v, w), the double bubble is completely determined by the mean curvatures 
of the outer caps. Thus, under-approximating the curvatures of these caps (by an 
analogous method to the single bubble case) corresponds to a double bubble with 
larger volumes and hence more surface area. 

We use Mathematica's built-in Set Accuracy function, which assures accuracy on 
25 binary digits right of the decimal point in our important computations. When we 
calculate curvature, we check that the calculated associated volume VolSphere[k] 
is less than v — 2"^'*, which guarantees that the associated volume is less that v. 
Hence the associated area will be less than A. Finally, we subtract 2~^* from our 
calculated area to be sure that we have a lower bound. 
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(IS^ 1/10,918^1/20) 




V 

Figure 6. The triangular domain of the {v,w) where numerical 
analysis is implemented. 



We check that the calculated associated volumes VolBubV[kl, k2] and 
VolBubW[kl, k2] exceed v and w by 2~^^ = 2 • 2~^'*', because sometimes there are 
two numbers summed in the calculation. Hence the associated area will be greater 
than A. We add 3 • 2^'^'* to the calculated area to be sure. 



For Claim 15.211 the proof-function is restricted to showing that the Hutchings 
function is positive only on rectangles that include the line v = .85w (See Figure 

Remark 5.2. In S^, the proof-function determines the size of each polygon via 
a recursive algorithm. Both proof-functions (Claims 15.51 and fS.Gp complete in well 
under six hours each depending on computer speed. In H^, the proof-functions 
(Claims l5. 10115. 20|) take well over 150 hours to complete. The function is broken 
into pieces in order to maximize efficiency by choosing the rectangle size to be 
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Figure 7. The steps needed to show that g > h on a triangular 
domain with hypotenuse along the line w = u. 
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Figure 8. Showing that the Hutchings function is positive on 
the rectangles covering the shaded region shows that the function 
is positive on the shaded region. 



relatively large on a given region. Finally, the claim dealing with the ray v = .85ui 
(Claim [5.2ip completes in under three hours. 
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5.2. Program implementation. This section discusses the code in the appendix. 
The computer proof in S'^ uses two proof functions. The first is called ProofFunc- 
tionTriangle. The other is called ProofPunctionRectangle. The following claims can 
be proved by examining and running the code in the Appendix (Section [7]). The 
proof-functions returns a 1 if the Hutchings function is positive on the rectangle or 
triangle it is called upon. 
Implementation in S^. 

Claim 5.3. The computed function A(v-, error-) returns a lower bound on 
the area of a sphere in with given volume v. 

Claim 5.4. The function A (v-, W-, VError-, WError-, changeToV-, change- 
ToW-) gives an upper bound on A(v,w). 

Claim 5.5. If ProofPunctionRectangle returns a 1, then the Hutchings Punc- 
tion is positive on this region. 

Claim 5.6. If ProofPunctionTriangle returns a 1, then the Hutchings Function 
is positive on this region. 

Claim 5.7. ProofPunctionRectangle[lhsF2, rhsF2, {VolOfS3/10, VolOfS3/10}, 
{VolOfS3/3, VolOfS3/3}, lhsP2[VolOfS3/10, 

VolOfS3/10], rhsP2[VolOfS3/10, VolOfS3/10], lhsP2[VolOfS3/3, VolOfS3/3], 
rhsP2[VolOfS3/3, VolOfS3/3j, temp] returns a 1 where 

lhsF2 =2* AreaSphereGivenV olume[v /2, V Error /2] + 
AreaSphereGivenV olume{w , W Error] + 
AreaSphereGivenV olume[v + w, V Error] 

and 

rhsP2 = 2* A[v, w, V Error, W Error, ChangelnV, ChangelnW]. 

Claim 5.8. Proof PunctionTriangle[lhsP2, rhsP2, VolOfS3/10, 
VolOfS3/3, VolOfS3/3, 9*VolOfS3/20, lhsP2[VolOfS3/10, 
VolOfS3/3], rhsP2[VolOfS3/10, VolOfS3/3], lhsP2[VolOfS3/10, 
9*VolOfS3/20], rhsF2[VolOfS3/10, 9*VolOfS3/20], 
lhsF2[VolOfS3/3, VolOfS3/3], rhsF2[VolOfS3/3, VolOfS3/3], 
temp] returns a 1 where 

lhsP2 =2* AreaSphereGivenV olume[v /2, V Error /2] + 
AreaSphereGivenVolume[w, W Error] 4- 
AreaSphereGivenV olume[v + w, V Error] 

and 

rhsP2 = 2* A[v, w, VError, W Error, ChangelnV, ChangelnW]. 
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Implementation in H^. The following claims can be proved by examining and run- 
ning the code in the appendix. The program stops if it finds a rectangle where 
g < h. Therefore, if the proof-function completes, the Hutchings function is posi- 
tive on the region it was called on (Claim [5^ . 

Claim 5.9. If ArrayFillingProof completes then the Hutchings function is pos- 
itive on the region it was called on. 

Claim 5.10. ArrayFillingProof [.002329, .00274, -01, .00001, 
.00001,11.46, 10.95, .9999, .9995] completes. 

Claim 5.11. ArrayFillingProof [.0085, .01, .1, .00005, .00005, 7.475, 7.15, 
.9999, .9995] completes. 

Claim 5.12. ArrayFillingProof [.085, .1, 1., .0005, .0005, 3.56, 3.415, .9999, 
.9995] completes. 

Claim 5.13. ArrayFillingProof [85, l.,15., .005, .005, 1.849, 1.787, .9999, 
.9995] completes. 

Claim 5.14. ArrayFillingProof [12.75, 15. ,25., .01, .01,1.15204, 1.1352, .9999, 
.9995] completes. 

Claim 5.15. ArrayFillingProof [21. 25, 25. ,45., .02, .02, 1.1027, 1.09054, -9999, 
.9995] completes. 

Claim 5.16. ArrayFillingProof [38. 25, 45. ,65., .02, .02, 1.0637, 1.05562, .9999, 
.9995] completes. 

Claim 5.17. ArrayFillingProof [55. 25, 65., 85., .02, .02, 1.04658, 1.040469, 
.9999, .9995] completes. 

Claim 5.18. ArrayFillingProof [72.25, 85., 110., .015, .015, 
1.03684, 1.031905, .9999, .9995] completes. 

Claim 5.19. ArrayFillingProof [93. 5, 110., 130., .015, .015, 
1.02927, 1.025276, .9999, .9995] completes. 

Claim 5.20. ArrayFillingProof [110.5, 130., 150., .015, .015, 
1.025161, 1.021693, .9999, .9995] completes. 

Finally if we adjust the ArrayFillingProof, so that it only shows that the Hutch- 
ings function is positive along the line v = .85w, 

Claim 5.21. ArrayFillingProof[127.5, 150., 300., .01, .01, 1.022077, 1.019009, 
.9999, .9995] completes. 

5.3. Main propositions. The following propositions consolidate the claims that 
the computer code runs showing the the Hutchings function is positive on the 
desired domains. Proposition [5211 uses a modified version of the code that checks 
only that the Hutchings function is positive on a small domain that includes the 
line V = .85ui. 

Proposition 5.22. If an area-minimizing double bubble in encloses volumes 
.1 < V < min{iL, 1 — 2w}, then the Hutchings function Fs3{v, w) > 0. 

Proof. This follows from implementing the code in the appendix (Claims [5.7|l5.8p , 
and Claims O and [ESI □ 
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Proposition 5.23. If an area-minimizing double bubble in tl^ encloses vol- 
umes v,w such that .002743 < v < ui < 150 and v > .85w then the Hutchings 
function Fii3{v,w) > 0. 

Proof This follows from Claim [5J] and Claims [ElOllEaOl □ 

Proposition 5.24. The Hutchings function Fiis{v,w) > 0, for volume pairs 
{v,w), where .85v = w and w £ [150,300]. 

Proof This follows from Claim [EH and Claim [5211 □ 
6. Double bubbles in and 

Theorems 16.31 and 16.41 provide our main results on double bubbles in and 
H^. They depend on Propositions 16.11 and 16.21 of Cotton and Freeman [CF| , af- 
ter Hutchings, Morgan, Ritore, and Ros [HMRR] , which reduce the proofs to our 
connectivity results of sections 3-5 via the Hutchings Function. 

Proposition 6.1. An area-minimizing double bubble in for which both en- 
closed regions and the exterior are connected must be standard. 

Proof. See Cotton and Freeman [CF[ Proposition 7.3]. □ 

Proposition 6.2. An area-minimizing double bubble in for which both 
enclosed regions are connected must be standard. 

Proof. See Cotton and Freeman [CF| Proposition 7.7]. □ 

Theorem 6.3 (Area-minimizing double bubbles in S^). An area-minimizing dou- 
ble bubble in is standard if the volume of each of the three regions is larger 
than 10% of the total volume of S^. 

Proof. By Proposition 16. 1| it suffices to show that all regions are connected. By 
Lemma [2. 14) we know that the region of volume v is connected whenever v > 2w 
01 V > 2u (Figure [9| a). This reduces the problem to showing that the region 
of volume v is connected on the pentagonal domain of Figure [9| e defined by the 
points: iv,in,u) = (.1, .1, .8), (.1, .8, .1), (.2, .7, .1), (.5, .25, .25), and (.2,.!,. 7). By 
Proposition 12. 3| it suffices to show that the Hutchings Function F{v,w) is pos- 
itive on this domain. When .1 < v < vj and v < 1 — 2vj, computer analysis 
(Proposition 15.221 see Figure [SJs) tells us that the Hutchings Function Fs3{v,w) 
is positive. In particular, for v < ^, Fs3{v,v) is positive. When ^ > v > w and 
V < 2w, by S-balancing (Proposition 12.131 see Figure[9l:) Fs3{v,w) is positive on 
the quadrilateral domain with endpoints: {v,w,u) = (.1, .1, .8), (|, |, |), (|, |, |), 
and (.2, .1, .7). Thus F{v,w) > on the pentagonal domain gained by S-balancing 
in addition to the region where the computer analysis was done. Hence by permu- 
tation (Proposition 12.151 see Figure [9ll), F{v,w) > on the pentagonal domain 
with endpoints: [v, w, u) = (.2, .7, .1), (|, i, |), (1, i, i), (.1, .45, .45) and (.1, .8, .1). 
In particular, if < .45, F{v,v) is positive. By S-balancing (Proposition 12. 13"1 see 
Figure[9t) F{v, w) is positive whenever in < v < 2w. Thus, F{v, w) is positive, for 
on the pentagonal domain {v,in,u) = (.1, .1, .8), (.1, .8, .1), (.2, .7, .1), (.5, .25, .25), 
and (.2,. 1,. 7). □ 

Theorem 6.4 (Area-minimizing double bubbles in H^). An area-minimizing 
double bubble in is standard if the volume of the smaller region is at least 
85% of the volume of the larger region. 





W V 

e) Balancing 
New Region; (1/3, 1/3, 1/3) 
(4/9, 1.3, 2/9), (.5, .25, .25), (4/9, 2/9, 1/3) 



Figure 9. An illustration of the steps taken to cover the whole 
domain of volumes on S^. 
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lim F = v=w 




V 



Figure 10. The basic steps of the proof that F is positive for 
V > .85w and w > .85v in H^. 



Proof. By Proposition 16.21 it suffices to show that the two enclosed regions are 
connected. By Proposition 12.3) it suffices to show that the Hutchings Function 
F{v,w) is positive for v > .85w (see Figure [T0|) . Computer analysis (Propositions 
[5?i3l [5T24|l shows that when .002743 < < w < 150 and v > .85w, the Hutchings 
Function Fii3{v, w) is positive. It also shows that the Hutchings function is positive 
along the line segment v = .85ii;,w £ [150,300]. Since the partial derivative of 
the Hutchings function with respect to v is positive for w > 150 (Lemma I4.22p . 
F{v,w) > for w e [150,300]. Asymptotic analysis (Theorems [3^ and l423l) shows 
that Fii3{v, w) is positive when v < w < .002743 or v/ .85 > ui > 300. In particular, 
Fhs{v,v) is positive. Hence when v >w, Proposition l2. 131 shows that Fii3{v,w) is 
positive. □ 
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Remark 6.5. In hyperbolic quotient spaces, Theorem 16.41 and (Cj Remark 4.4] 
prove that the standard double bubble is area minimizing for some small volume 
pairs where the smaller volume is at least 85 percent that of the larger. 

Future work. It follows from Theorem [631 and balancing (Lemma l2.14| Proposition 
I2.13P that the largest region in is always connected. Connectivity of the largest 
region is the key hypothesis for the even more general instability argument used in 
R'* [R, Section 8]; future work might produce an analog of this result for hyperbolic 
space. 

7. Appendix: program code 

7.1. code. 



(*The Volume of the 3-sphere with constant curvature 1*) 
VolDf S3-Pi» (2»Pi-Sin [2»Pi] ) ; 



(♦These functions give the relationships of illustrated in Figure 2.1*) 
ThetaFN [rl_ , r2.] : -ArcTan [Sqrt [3] * (Tan [r2] -Tan [rl] ) / (Tan [rl] +Tan [r2] ) ] ; 

xFN [rl_ ,r2_: : =ArcTan [Tan [rl] * (Sqrt [3] /2*Cos [ThetaFN [rl ,r2: ] +1/2* 

Sin [ThetaFN [rl , r2] ] ) ] ; 

rSFN [rl. , r2.] : -ArcCot [Cot [rl] -Cot [r2] ] ; 
PhilFN[rl.,r2J : -ArcSin[Sin[xFN [rl , r2] ] /Sin [rl] ] ; 

PsilFN[rl.,r2_] : -Module [{Phil} ,Phil-PhilFN[rl , r2] ;lf [(Tan[r2] -Tan[rl] )/( 
Tan[rl] +Tan[r2] ) -l/3<0 , Pi-Phil , Phil] ] 

Phi2FN[rl_,r2_] : -ArcSin[Sin[xFN [rl , r2] ] /Sin [r2] ] ; 

PhiSFN [rl_ , r2_] : =ArcSin[Sin[xFN [rl , r2] ] /Sin [r3FN [rl , r2] ] ] ; 



(♦Volume of a spherical cap in S3 of radius r and subtended by an angle 
of 2 Phi *) 

VolCap[r.,Phi.] : -If [r\ [Equal] Pi/2 , (* 

then*) ,(*else*)Pi (r-ArcTan [Cos [Phi] *Tan [r] ] -Cos [r] ♦Sin [r] +Cos [ 
Phi]*Cos[r]*Sin[r])] ; 

(*Area of a spherical cap in S3 of radius r and subtended by an angle 
of 2 Phi *) 

AreaCap[r.,Phi.] :-2Pi» (Sin [r] "2) * (1+Cos [Phi] ) ; 

(*Volume of a sphere of radius r in S3*) 
VolSphere [r_] : =Pi* (2*r-Sin [2*r] ) ; 

(*Area of a sphere of radius r in S3*) 
AreaSphereGivenRadius [r_] : =4*Pi* (Sin[r] ) ~2 

(*Standard double bubbles are compossed of two regions one enclosing the a 
smaller volume and a larger volume . The exterior cap of the 

smaller-volume-region has a radius of rl. The 

exterior cap of the larger-volume-region has radius r2.*) 

(*This function gives the volume of the smaller region of a standard double bubble 
determined by radii rl, r2*) 

VolDneSDB[rl_,r2J : -Module [{Psil , r3 ,Phi3> ,Psil-PsilFN [rl , r2] ; 
r3-r3FN[rl,r2] ; 
Phi3-Phi3FN[rl,r2] ; 

(*return*) VolCap [rl ,Psil] +VolCap [r3 ,Phi3] ] 

(*This function gives the volume of the larger region of a standard double bubble 
determined by radii rl, r2*) 

VolTuoSDB[rl_,r2.] : -Module [-CPhi2, r3 ,Phi3> ,Phi2-Phi2FN [rl , r2] ; 
r3-r3FN[rl,r2] ; 
Phi3-Phi3FN[rl,r2] ; 

(*retum*) VolCap [r2 , Pi-Phi2] -VolCap [r3 , Phi3] ] 

(*When both regions of the double bubble are the same it becomes advantageous 

to consider the one parameter function, that give the volume of one region given a 

radius r.*) 
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Equal VolGivenRadius [r_] : =Pi/2* ( 

2r-Siii[2r] ) * (1+ (Sqrt [2] *Cos [r] ) / CSqrt [7+Cos [2r] ] ) ) +Pi (ArcTan [Sqrt [2] Sin [r] \ 
/ CSqrt [7+Cos [2r] ] )] - (Sqrt [2] *r*Cos [r] ) / (Sqrt C7+Cos [2r] ] ) ) ; 



(♦Returns the sum of the areas of the three spherical caps that make up the standard 
double bubble corresponding to the two given radii*) 
AreaSDBGivenRadii[rl_,r2_] :=Module[{Psil,Phi2,Phi3,r3},Psil=PsilFN[rl,r2] ; 

Phi2=Phi2FN[rl,r2] ; 

Phi3=Phi3FN[rl,r2] ; 

r3=r3FN[rl,r2] ; 

(♦retum*) AreaCap [rl , Pi-Psil] +AreaCap [r2 ,Phi2] +AreaCap [r3 ,Pi-Phi3] ] 

(^Returns a radius of a sphere in S3 has less area than the sphere of volume v. 
If V <= VolDfS3/2, then this function returns a radius of a sphere with volume in 

the interval [v-error, v] . If v > VoiafS3/2, then 
this function returns a radius of a sphere with volume in the interval [v, v+error] *) 
RadiusSphere [v_ , error_] : =Module [-CminRadius , mazRadius , currentRadius , 
counter} , 

counter =0; 
minRadius=0 ; 
maxRadius = Pi; 
currentRadius = Pi/2; 

If [v VolQfS3/2, 

(♦This while loop checks to see if the volume is in the interval [v-error, v] if not 
then it checks to see if it the volume is two high or low. The While loop 
improves the estimates minRadius and maxRadius, While the loop is running 
minRadius corresponds to sphere of volume less than v-error and maxRadius 
corresponds to a sphere of volume bigger than v. 

Finally, the sphere of radius currentRadius corresponds to a sphere of volume 
in [v-error, v] and the loop stops.*) 

While [VolSphere [currentRadius] > v II (VolSphere [currentRadius] < v \ 

-error) , 

If [VolSphere [currentRadius] > v, 
iiiaxRadius= currentRadius ; 

currentRadius = ( currentRadius + minRadius ) /2 ; 

]; 

If [VolSphere [currentRadius] < v-error , 
minRadius= currentRadius ; 

currentRadius = ( currentRadius + maxRadius ) /2 ; 

]; 

count er++; 
]. 

(♦else if V > volofs3/2^) 

(♦This while loop checks to see if the volume is in the interval [v, v+error] if not 
then it checks to see if it the volume is two high or low. The While loop 
improves the estimate minRadius and maxRadius. While the loop is running 
minRadius corresponds to sphere of volume less than v and maixRadius 
corresponds to a sphere of volume bigger than v+error. Finally, the sphere of 
radius currentRadius corresponds to a sphere of volume in [v, v+error] 
and the loop stops. ♦) 

While [VolSphere [currentRadius] < v || (VolSphere [currentRadius] > v \ 

+error) , 

If [VolSphere [currentRadius] < v, 
minRadius= currentRadius ; 

currentRadius = (currentRadius + maxRadius) /2 ; 

]; 

If [VolSphere [currentRadius] > v + error, 
maxRadius= currentRadius ; 

currentRadius = (currentRadius + minRadius) /2; 
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]; 

counter++ ; 
]; 

]; 

(♦Print ["It took " , counter, "steps . "] ; 

Print ["Radius is " , N [currentRadius] ] ; *) 
currentRadius 

h 



(♦Gives an overestimate for the sphere of volume v*) 

AreaSphereGivenVolume [v_, error _] :=AreaSphereGivenKadius [RadiusSphere [v, \ 
error]] ; 

(* This function returns a pair of radii corresponding to a standard double bubble 
with volumes bigger than v,w by making very small adjustments to the radii of the 
exterior caps. This function will only work for v<= w < u. *) 

RadiiSDB [v_ , w_ , vError_ , wError_ , vChangingFactor_ , 

wChangingFactor_] : =Module [{VRadius , WRadius , madeChangeToV , \ 

madeChangeToW , ad j ust edVErr or , 

adjust edWError , vAdjustment , wAdjustment .volOne , volTwo, counter} , 

VRadius = Pi/4; 
WRadius = Pi/4; 
count er=0; 



vAdjustment = vCheingingFactor ; 
wAd j ustment =wChangingFact or ; 

(♦These variables ma^e sure the v, w we approximate are in the region 
we where A(v,w) increases in v and w*) 
adjustedVError = vError; 
adjust edWError = wError; 

(♦assurance that we are below line w =u and above line v=w^) 
If [v+2w >= VolDfS3 II v >= w, 

Print["ERRDR RadiiSDB: v and w are not in increasingRegion! "] ; 

Print ["v+ 2w = " ,N[v+2w]]; 

Print ["Vol of S3 = N[Vol0fS3]]; 

Print ["v =" ,v] ; 

Print ["w =" ,w] ; 



While [v+adjustedVError + 2(w+adjustedWError) > Vol0fS3 II \ 
v+adjustedVError > w+adjustedWError, 

If [v+adjustedVError > Vol0fS3 - 2+w, 

adjustedVError= adj stedVError/2 ; 

]; 

If [2(w+adjustedWError) > Vol0fS3 - (v + adjustedVError), 

adj ustedWError= adj ustedWError/2 ; 
adjustedVError= adjustedVError/Sqrt [2] ; 

]; 

If [v+adjustedVError > w+adjustedWError, 
adjustedVError = adjustedVError/2; 

]; 

]; 

While [(VolOneSDB [VRadius, WRadius] <v || (VolOneSDB [VRadius , 
WRadius] > v+ 

adjustedVError) I I VolTwoSDB [VRadius , 

WRadius] < w I I (VolTwoSDB [VRadius , WRadius] > \ 
w+adjustedWError)) , 



(♦These two variable prevent infinite while loop from occxiring*) 
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madeChangeToV = False; 
madeChangeToW = False; 

If tVolOneSDBEVRadius, WRadius]< v, 

VRadius = VRadius*vAdjustment; 
iiiadeChaiigeToV= ! madeChangeToV ; 

]; 

(♦This set of if statements adjusts one or both of the radii. If it can 
not change one of the radii then it refines the amount that the radii 
are adjusted by,*) 

If [VolOneSDB [VRadius, WRadius]> v+adjustedVError, 

VRadius = VRadius/vAdjustment ; 
madeChangeToV= ! madeChangeToV ; 

]; 



If [VolTwoSDB [VRadius , WRadius]< w, 



MRadius = WRadius*wAdjustment ; 
madeChangeToW= ! madeChangeToW ; 

]; 

If [VolTwoSDB [VRadius , WRadius] > w+adjustedWError , 

WRadius = WRadius/wAdjustment ; 
iiiadeChangeToW= ! madeChangeToW ; 

]; 

If [! madeChangeToV, 

vAdjustment=CvAdjustment+l) /2 ; 

]; 

If [! madeChangeToW, 

wAdjustment=CwAdjustment+l) /2 ; 

]; 

(♦Print ["Volume Vat ", N [VolGneSDB [VRadius , WRadius]]]; 
Print ["Volume W at ", N [VolTwoSDB [VRadius , WRadius]]]; 

Print ["While Loop test 

is ", (VolQneSDB [VRadius, WRadius] <v || (VolOneSDB [VRadius , \ 
WRadius] > v+adjustedVError) I I VolTwoSDB [VRadius , WRadius] < 
w II (VolTwoSDB [VRadius, WRadius] > 
w+adjustedWError))] ;*) 

]; 

]; 

(♦Print["It took this many steps: ", counter] ;♦) 

{VRadius, WRadius} 

]; 

(♦By symmetry the standard double bubble enclosing volumes (v,w,u) has the same 
area of the standard double bubble enclosing volumes (u,w,v). So, when w=u we can 
get an overestimate for the area of the stajidard double bubble enclosing volumes 
(VjUjU) by geting an overestime for the area of the standard double bubble 
enclosing volumes (u,w,u) . This function returns a pair of radii that correspond to a 
double bubble where the larger volume is equal to the exterior and 
both are bigger than or equal to w*) 
RadiiWhenWEqualsU[w_, 
error _] : =Module [{radius , minRadius , maxRadius , counter , adjustedError} , 
minRadius =0 ; 
maxRadius=3Pi/2 ; 

radius = (minRadius+maxRadius/2) ; 
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adjust edError=error; 
count er=0; 

If[3w< Vol0fS3, 

Print ["W is too small,"], 
While[3(w-adjustedError)< VolQfS3, 

adjustedError = adjustedError/2; 

]; 

C*Print ["Target Range - C", N [u-error] , " , " , N[w],")" ];*) 

While [EqualVolGivenRadius [radius] >w I I 
EqualVolGivenRadius [radius] <w-adjustedError , 
C*counter++; 

Print ["counter = ", counter] ;*) 



If [EqualVolGivenRadius [radius] >w, 

maxRadius =radius ; 

radius = CiiiinRadius+radius)/2; 

]; 

If [EqualVolGivenRadius [radius] <w- adjustedError , 

minRadius =radius ; 

radius = (maxRadius+radius)/2; 

]; 

C*Print["Vol = " ,W[EqualVclGivenRadius [radius] ]]; *) 

] ; 

C*Print["It took this many steps: ", counter];*) 

]; 

{radius , radius} 
]; 



C*This piece of code should be used when v=w*) 
RadiiWhenVEqual sW [v_ , error_] : =Module [{radius , 
minRadius , maxRadius , counter , adjustedError} , 

minRadius =0 ; 

maxRadius=Pi ; 

radius = CminRadius+maxRadius/2) ; 
adjust edError=error; 
count er=0; 

(♦Print ["radius= " , radius] ;*) 
If[3v>= VolOfSa, 

Print ["V is too big."], 

While [3 (v+adjustedError)> VolDf S3 , 
adjustedError = adjustedError/2; 

]; 

(♦Print ["Target Range = C", N [v-error] , " , " , N[v],")" ];*) 

While [EqualVolGivenRadius [ 
radius] <v I I EqualVolGivenRadius [radius] >v+adjustedError , 



(♦counter++; 

Print ["counter = ", counter] ;♦) 



If [EqualVolGivenRadius [radius] >v+adj ustedError , 
maxRadius =radius ; 
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radius = (minRadius+radius) /2 ; 

]; 

If [Equal VolGivenRadi US [radius] <v, 

minRadius =radius ; 

radius = CmaxRadius+radius)/2; 

]; 

C*Print["Vol = " ,N[EqualVolGiveiiRadius [radius] ]]; *) 
]; 

C*Print["It took this many steps: ", counter];*) 

]; 

{radius , radius} 



(♦This function returns an overestimate for the area of a 
standard double bubble of volumes v, w.*) 

A[v_,w_, vError_, 

wError_ , vChaiigingFactor_ , \ 
wChangingFactor_] : =Module [{radii ,u, rl , r2 , area} , 

u=Vol0fS3-v-w; 

(♦when all three volumes are equal*) 
If [v==w && w==u, 

area = 6*Pi; 



(♦one the line v=w^) 
If [v==w, 

radii = RadiiWhenVEqualsW[v, vError] ; 

area = AreaSDBGivenRadii [radii [ [1] ] , radii[[2]]]; 



(♦on the line w=u^) 
If [w==u, 

radii = RadiiWhenWEqualsU[w, wError] ; 

area = AreaSDBGivenRadii [radii [ [1] ] , radii [[2]]]; 



(♦If all volumes are distinct*) 

radi i = Radi iSDB [v , w , vError , wError , vChangingFact or , 
wChangingFactor] ; 
area = AreaSDBGivenRadii [radii [ [1] ] , radii [[2]]]; 

]; 

]; 

]; 

area 
]; 

(♦This function is described in Section 4 of the paper. It shows that either 
the Hutchings function is on a rectangular domain or it breaks up the 
rectangle into 4 smaller rectangles and checks the code again. The function 
completes and returns a 1 if the function is positive on the rectangular domain 
and if the approximated function is negative on any part of the domain^) 
Proof FunctionRectangle [lhs_,rhs_,pl_,p3_,leftll_, right ll_,left33_, 

right 33_ , depth.] : =Module [{decision , xl , y 1 , x2 , y2 , xS , 

y3, left 12, rightl2, left 13, right 13, 

lef t21 ,right21 , lef t22 , right22 , lef t23, right23 , lef t31 , 
right31 , lef t32 , right 32 , dl , d2 , d3 , d4} , depth=l ; 



xl=pl[[l]]; 
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yl=pl[[2]] ; 

x3-p3[[l]] ; 
y3=p3[[2]]; 

x2=Cxl+x3)/2; 

y2=Cyl+y3)/2; 

leftl2=lhs [xl,y2] ;rightl2=rlis [xl,y2] ; 
Ieftl3=lhs[xl,y3] ;rightl3=rhs [xl,y3] ; 
Ieft21=lhs[x2,yl] ;right21=rhs [x2,yl] ; 
lef t22=lhs [x2 , y2] ; right22=rhs [x2 , y2] ; 
lGft23=lhs Cx2,y3] ; right23-rhs [x2,y3] ; 
lef t31=lhs [x3 , yl] ; right31=rhs [x3 , yl] ; 
Ieft32=lhs[x3,y2] ;right32=rhs [x3,y2] ; 

(*if xl > y3 then balancing 

covers the region and we don't need to check it.*) 
If[xl > y3, 

decision =1, 

(♦otherwise we do*) 

If [left 1 l>right 1 Ifc&lef t33>right33 , 



(*then*) 

If [Min[leftll,leftl3,left31,left33]>right33, C*then*)Print ["Points \ 
",pl," and ",p3," — Direct hit!"]; 
decision=l, 

(*else*) Print ["Points " ,pl, " and " , 
p3," — Splitting into four."]; 

If [Proof FunctionRect angle [Ihs,rhs,-Cxl,yl},{x2,y2}, lef til, right 11, 
lef t22,right22,dl]\ [Equal] 1&& 

Proof FunctionRectangle [Ihs , rhs , {x2 ,yl} , {x3 ,y2} , lef t21 , right21 , \ 
Ieft32,right32,d2]\ [Equal] 1&& 

Proof FunctionRectangle [ 

Ihs, rhs, {xl,y2>,-[x2,y3}, lef tl2, right 12, lef t23,right23, 

d3] \ [Equal] l&fcPr oof FunctionRectangle [Ihs , rhs , ■[x2 , y2}, { 
x3, y3} , lef t22, right22 , lef t33 ,right33 ,d4] \ [Equal] 1 , 

C*then*) 

Print ["Points ",pl," and 

",p3," — Hit after checking four! Depth: ", Max [depth, 

dl+1 , d2+l , d3+l , d4+l] ] ; 
decision=l , 

C*else*)Print ["Proof function failed"] ; 

decision=0 ; 
C*endif*)] ; 

depth=Max [depth, dl+1 , d2+l , d3+l , d4+l] ; 

(*endif*)] .Print ["Oh no for ",pl," and ",p3,"!"]; 
Print [N [lef til]] ; 
Print [N[rightll]] ; 
Print [N [lef t33]] ; 
Print [N [right33] ] ; 
decision=0 ; 
]; 

(*endif *) 
]; 

C *r et urn* ) dec! s i on 
]; 

C*This function is described in Section 4 of the paper. It shows that either 
the Hutchings function is on a triangular domain or it breaks up the 

triangle into 2 smaller triangles and a rectangle and checks the code again. The function 
completes when it has broken the domain into small enough domains to see that the function 
is postive and returns a 1 or if the approximated function is negative on any part of the 
domain the function returns a 0.*) 

ProofFunctionTriangle [Ihs_,rhs_,xl_,yl_,x3_,y3_, lef tll_, right 11_, lef tl3_,\ 
right 13_, lef t31_,right31_, depth.] :=Module [{decision, x2,y2, left 12, right 12, \ 
lef t21 ,right21 , lef t22 , right22 , dl , d2 ,d3} ,depth=l ; 

Print ["starting up the process"]; 



x2=Cxl+x3)/2; 
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y2=(yl+y3)/2; 
Print ["whoa now!"]; 

If [Ieftll>rightll&&leftl3>rightl3&&left31>right31, C*then*) 
H[Min[leftll, 

lef tl3,lef t31] >right31 , C*tli8n*)Print ["Points " ,yl}, " , " ,-[xl , 
y3},", and ",{x3,yl}," — Direct triangle hit!"]; 
decision=l, 

C*else*)Print ["Points ",{xl,yl},", " ,-[xl ,y3}, " , and ",{x3,yl}," — 
Splitting into three."]; 

Ieftl2-lhs [xl,y2] ;rightl2-rhs [xl,y2] ; 
Ieft21"lhs[x2,yl] ;right21"rhs [x2,yl] ; 
lef t22=lhs [x2 , y2] ; right22=rhs [x2 , y2] ; 

If [Proof FunctionRect angle [Ihs , rhs , {xl , yl} , -[x2 , y2> , left 1 1 , right 11 , \ 
lef t22,right22,dl]\ [Equal] l&&Proof FunctionTriangle [Ihs ,rhs ,xl ,y2 , x2 ,y3,lef tl2, 

rightl2,leftl3,rightl3,left22,right22,d2]\ [Equal] 
Proof FunctionTriangle [Ihs , rhs , x2 , yl , x3 , y2 , lef t21 , 

right21 , lef t22 , right22 , lef t 31 , rightSl , d3] \ [Equal] 1 , C*then*) 
Print ["Points ",{xl,yl},", ",{xl,y3>,", and ",-[x3,yl>," — 
Hit after checking three! Depth: " ,Max[depth,dl+l ,d2+l ,d3+l] ] ; 
decision=l, C*else*)Print [ 

"Points ",-(xl,yl},", " ,-(xl,y3}, " , and ",-[x3,yl}," — Failed 
after checking three!"]; 
decision=0; 
C*endif *)] ; 
depth=Max[depth,dl+l,d2+l,d3+l] ; 
(*endif *)] , C*else*)Print ["Oh no for " -(xl.yl},", 
",{xl,y3>,", and " ,{x3,yl}, " ! "] ; 
Print [N [lef til]] ; 
Print [»[rightll]] ; 
Print [B [lef tl3]] ; 
Print [N [rightlS] ] ; 
Print [N [lef t31]] ; 
Print [N[right31] ] ; 
decision=0; 
(*endif *)] ; 
(♦return*) decision] 



H3 Code 

C*error control in computation the default value of maching*) 
WorkingPrecision->MachinePrecision; 
C0MPUTE31E3U10R= 2--24; 
ACCURACY = 2B; 



LargestSingCurvature = 16.8; 



AreaSphere [k_Real] ;=SetAccuracy [C4*Pi)/ C-l+k"2)-C0MPUTERERR0R, ACCURACY] ; 

VolSphere [k.Real] : -Set Accuracy [Pi* (-2* 

ArcCoth [k] +Sinh [2*ArcCoth [k] ] ) , ACCURACY] ; 

ASClclose [kl_Real , k2_Real] : =SetAccuracy [ C2*Pi* Cl+Sqrt [- C Ckl"2* C-l+Cos [ Cl/6) * C\ 
Pi-6*ArcTan [ (Sqrt [3] * (kl-k2) ) / ( 

kl+k2)])]-2))/( 

kl-2-Cos[(l/6)*( 

Pi-6*ArcTan[C 
Sqrt [3]*Ckl- 

k2))/(kl+k2)])]-2))]))/(-l+kl-2). \ 

ACCURACY+1] ; 
ASClfar[kl_Real, 

k2.Real] ; =SetAccuracy [- ( (2*Pi» (-1+Sqrt [- ( (kl-2* (-1+Cos [ (1/6) * (Pi-6*\ 
ArcTan [ (Sqrt [3] * (kl-k2) ) / ( 

kl+k2)])]-2))/(kl-2-Cos[( 

l/6)*(Pi-6*ArcTan[( 

Sqrt[3]*( 

kl-k2))/(kl+k2)])] -2) )]))/(-!+ 
kl-2)), ACCURACY+1]; 

ASCl[kl_Real,k2_Real] : -Set Accuracy [If [ (kl-k2) / (kl+k2)-l/3<0, ASClclose [ 
kl,k2] ,ASClfar[kl,k2]] , ACCURACY+1] ; 

ASC2 [kl.Real ,k2_Real] : =SetAccuracy [ (2*Pi* ( 1+Sqrt [ (kl-2-k2-2*Cos [( 1/6) * (Pi-6* 
ArcTan[(Sqrt [3] * (kl-k2) ) / (kl+ 
k2)])]-2)/( 
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kl-2-Cos[(l/ 

6) * CPi-6*ArcTan [ (Sqrt [3] » (kl-k2) ) / (kl+ 

k2)])]-2)]))/(-l+k2-2), ACCURACY+1]; 

ASC3 [kl.Real ,k2_Real] : =SetAccuracy [ (2*Pi* (1-Sqrt [(kl-2- (kl-k2) -2*Cos [ ( 
l/6)*CPi-6* 

ArcTaii[(Sqrt[3]*C 
kl-k2))/(kl+ 

k2)])r2)/(kl-2- 
Cob[(1/6)*( 
Pi-6*ArcTan [ C 
Sqrt[3]*(kl- 

k2))/(kl+k2)])]-2)]))/(-l+(kl-k2)-2), \ 

ACCURACY+1] ; 

AreaDblBubble [kl.Real ,k2_Eeal] : -Set Accuracy [ASCI [kl ,k2] +ASC2[kl , 
k2]+ASC3[kl,k2]+3*C0MPUTERERRaR, ACCURACY] ; 

VolCaplclose [kl.Real , k2.Real] :=-((!/ (-l+kl-2) ) * CPi* ( C-l+kl-2) * 

ArcCotli[kl] + C-l+kl-2)»ArcTanh[Sqrt [-( (kl-2»C-l+Cos [(1/6)* (Pi-6*ArcTaii[(\ 
Sqrt [3] * (kl-k2) ) / Ckl+k2) ] ) ] - 2) ) / (kl -2-Cos [Cl/6) « (Pi-6*ArcTaii [ (Sqrt [3] * (kl-k2) ) \ 
/ Ckl+k2) ] ) ] -2) ) ] /kl] -kl* ( 1+Sqrt [- ( ( 

kl-2* (-1+Cos [Cl/6) * (Pi-6*ArcTan [(Sqrt [\ 
3] * (kl-k2) ) / (kl+k2)] ) ] -2) ) / (kl-2-Cos [ (1/6) * (Pi-6*ArcTaii [ ( 

Sqrt [3] * (kl-k2) ) / (kl+k2) ] ) ] -2) ) ] ) ) ) ) ; 

VolCaplfar [kl.Real, k2_Real] : "Pi* (-ArcCoth[kl]+{l/ (-l+kl-2) ) *(kl+(-l+kl-2) *\ 
ArcTanh [Sqrt [- ( (kl-2* (-1+Cos [(1/6) * CPi-6*ArcTaii [(Sqrt [ 
3]*(kl-k2))/( 

kl+k2)])]-2))/(kl-2-Cos[(l/6)*(Pi-6*ArcTaii[(Sqrt[ 
3]*(kl- 

k2))/(kl+ 

k2)])]-2))]/ 
kl]-kl» 

Sqrt [-C (kl-2* (-1+Coa [(1/6)* (Pi-6*\ 
ArcTan [ (Sqrt [3] * (kl-k2) ) / (kl+k2)] )] -2) ) / (kl-2-CoB [ ( 

1/6) * (Pi-6*ArcTan [ (Sqrt [3] * (kl-k2) ) / (kl+ 

k2)])]-2))])); 

VolCapl [kl.Real , k2.Real] : "If [(kl-k2) / (kl+k2) -l/3<0 , 
VolCaplclose [kl,k2] , VolCaplfar [kl, k2]] ; 

VolCap2 [kl.Real , k2.Real] : "Pi* (- (k2/ (l-k2-2) ) -ArcCoth [k2] - 

ArcTanh [Sqrt [(kl-2-k2-2»Cos [ Cl/6) * (Pi-6* ArcTan [(Sqrt [3] *( 
kl-k2))/(kl+k2)])]-2)/(kl-2-Cos[( 
1/6)* (Pi-6* 

ArcTan [(Sqrt [3] *( 

kl-k2) ) / (kl+k2) ] ) ] -2) ] /k2] - (k2*Sqrt [ (kl-2-\ 
k2-2*Cos [ (1/6) * (Pi-6*ArcTan [ (Sqrt [3] * (kl-k2) ) / Ckl+k2)] )] -2) / (kl-2- 
Cos [ (1/6) * (Pi-6*ArcTan [ (Sqrt [ 

3] * (kl-k2) ) / (kl+k2)] )] -2)] ) / (l-k2-2) ) ; 

VolCapS [kl.Real, 

k2_Real] : =Re [Pi* (ArcTanh [Sqrt [ (kl-2- (kl-k2) -2«Cos [(1/6) * (Pi-6*ArcTan [ (\ 
Sqrt [3] * (kl-k2) ) / (kl+k2) ] )] -2) / (kl-2- 

Cos [ (1/6) * (Pi-6»ArcTan [ (Sqrt [ 

3]*(kl-k2))/(kl+k2)])]-2)]/(kl-k2)]+(l/(-\ 
l+kl-2-2*kl*k2+k2-2) ) * (kl-k2-(-l+kl-2-2* 
kl*k2+k2-2)* 
ArcCoth [ 

kl-k2]+(-kl+ 

k2)*Sqrt [(kl-2- (kl-k2)-2*Cos [(1/6)* (Pi-6*ArcTan[\ 
(Sqrt [3] »(kl-k2))/(kl+k2)] )] -2)/ (kl-2-Cos [(1/6) * CPi-6*ArcTan[C 
Sqrt [3]»Ckl-k2))/Ckl+k2)])]-2)]))]; 

VolBubl[kl.Real,k2.Real] :-If [kl < 1 II k2 < 1, 

Print ["VolBubl Error : kl or k2 too small"]; 
Print ["kl " " , kl , " k2 " " , k2 ] ; 
,If[kl \[Equal] k2,VolCapl[kl,k2] 

VolCapl [kl , k2] +VolCap3 [kl , k2] 
] 

]; 



VolBub2 [kl.Real, k2.Real] : "If [kl\ [Equal] k2, VolCap2 [kl ,k2] 

VolCap2 [kl ,k2] -VolCap3 [kl ,k2] 
]; 

VolBublA [kl.Real, k2.Real] :=If [kl>" k2, 
VolBubl [kl, k2] 
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C*Print ["VolBublA curvatures switched"];*) 

VolBub2[k2, kl] 
]; 

VolBub2A [kl_Real , k2_Real] :=H [kl>= k2, 
VolBub2[kl, k2] 

C*Print ["VolBub2A curvatures switched to ", VolBubl[kl, k2]] ;*) 

VolBubl[k2, kl] 

]; 

VolBubV[kl_Real, k2_Real] := SetAccuracy [VolBublA [kl, k2] , ACCURACY]; 



VolBubW[kl_Real, k2_Real] := SetAccuracy [VolBub2A[kl, k2] , 16]; 

C*gives an underappoxiamation for the curvature of a sphere in H3*) 
CurvatureEstimatorSingle[v_Real] := If [v > .0001, 
If [ V < .008, 
8.109, 

If[ v<.02, 

6. , 

If [ 

(*.02 V < 1*) 
v < 1, 

1.84871, 
If [ 

(•1 V < 10*) 

v<10, 

1.19394, 

If [ 

(»20 v < 30 *) 
V < 30, 
1.08108, 

If [ 

C*30 V < 50*) 
v<BO, 

1.05231, 

If [ 

(*60 V < 60*) 
v< 60, 

1.04466, 

If [ 

(»60 <- V < 70*) 
v< 70, 

1.03883, 

If [ 

C*70 V < 80*) 
v< 80, 

1.03437, 

If[ 

C*80 V < 90*) 
v< 90, 

1.03086, 

If [ 

(»90 v < 110*) 
v<110. 



1.0266, 
If[ 

(*110< V < 200*) 
V < 200, 
1.01468, 
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(♦else 550>v 110*) 

1.00B09 

] 

] 

] 

] 

] 

] 

] 

] 

] 

] 

] 

], 

C*else if V < .0001*) 
Print ["volume to small"]; 
C*return*) 
-1 

]; 

(♦This function uses interpolation to get asure that the curvature used in \ 
the approximation is less than the given v within 

an interval that is small. The constant term COMPUTERERRDR assures that \ 

the precision of the computer does not affect our lower bound*) 
Curvaturef romVolSingle [v_Real , error_Real] : ^Module [{ 

curvature, smallK = CurvatureEstimatorSingle [v] , bigK= \ 
LargestSingCurvature, counter =0}, 

curvature = smallK; 

(**) 

While [VolSphere [curvature] > v - COMPUTERERROR I I VolSphere [curvature] \ 
< v-error , 

curvature = CbigK+smallK)/2; 
counter = counter +1; 

C*this tightens the interval that curvatures can be in*) 
If [VolSphere [curvature] > v - COMPUTERERROR, 

smallK= curvature ; 

, If [VolSphere [curvature] < v - error, 
bigK = curvature; 

]; 

] ; 

] ; 

(♦Print C"It took " , counter, "steps. " ] ;♦) 
curvature 

]; 



(♦this function returns curvature pair that is at least one box size over \ 
then returns the array with the new curvature pair and the changes for each \ 
curvature though this seems to do the same thing as Curvaturesf romVolDouble [] \ 
this doesn't limit the size of the change^) 

MakeChangeOneBoxSizeDver [kl_Real , k2_Real , volume l_Real , BoxSize_Real , \ 
changel_Real] : = Module [{minVolCheinge , maxVolChange , ad jCurvaturel , \ 
adjCurvature2, ad j Change 1 , f inalArray , counter, slope , closest OverShot , \ 
close stUnderShot , lastCurvatureGuess}, 

minVolChange = . 8*BoxSize ; 
maxVolChange = 1.7*BoxSize; 

adjCurvaturel = changel^kl; 
adjCurvature2 = k2; 
ad j Change 1 = change 1 ; 
count er=0; 

closestOverShot= k2 ; 
closestUnderShot = kl; 

slope = (VolBublA [adjCurvaturel, adjCurvatur82] - \ 
volumel)/(adjCurvaturel - kl) ; 

adjCurvaturel = adjCurvaturel + BoxSize/ slope ; 

(♦Print ["curvatures, kl = ", adjCurvaturel, " k2 is ", k2] ;♦) 
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WhileC CVolBublA[adjCurvaturel,adjCurvatiire2]- volumeK minVolChangG I I \ 
VolBublA[adjCurvaturel,adjCurvature2]- volumel > maxVolChange) ftft \ 
adjCurvaturel > 1 M counter < 10, 

If [VolBublA [adjCurvaturel ,adjCiirvature2] - volumeK mlnVolChange , 

C*Print["hi test 1 is VolBublA[ 

adjCurvaturel, ad jCurvature2]- volumeK mlnVolChange] ; *) 
lastCurvatureGuess = adjCurvaturel; 

adj Curvaturel = (adj Curvaturel + closest OverShot ) /2 ; 
closestOverShot = lastCurvatureGuess ; 



If [VolBublA [adjCurvaturel , adj Curvature2] - volumel> mazVolChange , 

(*Print ["hi test 2 is " , VolBublA [adjCurvaturel , adj Curvature2] - 
volumel> maxVolChange] ; *) 

lastCurvatureGuess = adjCurvaturel; 

adjCurvaturel = (adjCurvaturel + closest UnderShot) /2; 

closestUnderShot = lastCurvatureGuess ; 
]; 

C*Print ["trial" , counter , " vol is \ 
" ,VolBub [adjCurvaturel , adj Curvature 2] ] ; *) 

counter^ counter+1 ; 
(♦Print ["counter = ", counter] ; 
Print ["voll =", 

VolBublA [adjCurvaturel , adj Curvature2] ] ; 
Print ["change = \ 
" , adjChangel] ; 

Print ["curvature " , adjCurvaturel ] ; 

*) 

]; 

If [adjCurvaturel <= 1 , 

Print ["MakeChangeQneBoxSizeOver error curvature too small"] ; 

]; 

C* 

Print["voll =" , VolBubl [adjCurvaturel, adjCurvature2] ] ; 
Print ["change = ", adjChangel] ; 
Print ["curvature adjCurvaturel] ; 
*) 

adjChangel = adjCurvaturel/kl ; 

{adjCurvaturel , adjChangel} 
]; 



(*this function returns curvature pair that is at least one box size over \ 
then returns the array with the new curvature pair and the changes for each \ 
curvature though this seems to do the same thing as CMrvaturesf romVolDouble [] \ 
this doesn' t limit the size of the change*) 

MakeChangeGneBoxSizeUp[kl_Real, k2_Real, volume2_Real , BoxSize_Real , \ 
change2_Real] : = Module [{minVol Change , maxVolChange , adjCurvaturel , \ 
adj Curvature2 , adjChangel , counter , slope , lastCurvatureGuess , \ 
closestUnderShot , closestOverShot} , 

minVolChange = l.*BoxSize; 
maxVolChange = 2.*BoxSize; 

adjCurvaturel = kl; 

ad j Curvature2 = k2*change2; 

adjChange2 = change2; 

close stUnderShot=k2; 

ClosestOverShot = 1.; 

(♦Might want to change this!*) 

counter=0; 

C*Print["hi"] ; 

Print ["test =" , (VolBub2 [adjCurvaturel , adj Curvature2] - volume2< \ 
mlnVolChange I I VolBub2[adjCurvaturel, 
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adjCurvature2] - voliime2 > maxVolChange) && 
adjCurvature2 > 1 ftft counter < 10] ; 

*) 

slope = C 

VolBub2A[adjCurvaturel, adjCurvature2] - voliime2)/CadjCurvature2 - k2) ; 

adjCurvature2 = adjCurvature2 + BoxSize/ slope ; 

While[ CVolBub2A[adjCurvaturel, 
adjCurvature2] - voliLme2< minVolChange I I \ 
VolBub2A [ad j Curvature 1 , adjCurvature2] - 

volume2 > mazVolChange) &ft ad j Curvatiire2 > 1 &ft counter < 10, 

counter= counter+1 ; 

(*Print ["counter = " , counter] ; 

Print["vol2 =", VolBub2[adjCurvaturel,adjCurvature2] ]; 

Print ["change = adjChaiige2] ; 

Print ["curvature " , adjCurvaturel ] ; 

*) 

If [VolBub2A[adjCurvaturel ,adjCurvature2] - volumeK minVolCliange , 
lastCurvatureGuess = adjCurvature2; 

adjCurvature2 = CadjCurvature2 + closestOverShot)/2; 

closestOverShot = lastCurvatureGuess; 
]; 



If [VolBub2A [adjCurvaturel , adj CurvatureS] - volume 1> raaxVolCliange , 
lastCurvatureGuess = adjCurvaturel; 

adjCurvature2 = CadjCurvature2 + closest UnderShot) /2 ; 
closestUnderShot = lastCurvatureGuess ; 



]; 

(♦Print ["counter = " , counter] ; *) 
If [adjCurvature2 <= 1, 

Print ["MakeChangeOneBoxSizeUp error curvature too small"]; 

]; 

C* 

Print["vol2 =" , VolBub2 [adjCurvaturel ,adjCurvature2] ]; 
Print ["change = " , adjChange2] ; 
Print ["curvature " , ad j Curvature2] ; 
*) 

adjChange2 = adjCurvature2/k2; 
■[adjCurvature2 , adjChange2} 



C*This function tells whether the point is in given rectangle with lower left \ 
corner xvalueBox, yValueBox and width boxWidth and heigh boxHeight*) 
IsPointinBox [xValuePoint_Real , yValuePoint_Real , xValueBox_Real , \ 
yValueBox_Real , boxWidth_Real , boxHeight_Real] : = 

(xValuePoint> xValueBox+ COMPUTERERROR && xValuePoint < (xValueBox + \ 
boxWidth) && yValuePoint > yValueBox + COMPUTERERROR ftft 
yValuePoint < (yValueBox + boxHeight)); 



(♦This function tells whether the point has a x value and a yalue bigger than \ 
a given x, and y value ♦) 

IsPointToUpperRight [xValuePoint_Real , yValuePoint _Real , xValueBox_Real , \ 
yValueBox^Real] : = 

(xValuePoint> xValueBox) && yValuePoint > yValueBox ; 

(♦This function is supposed manipulate a given curvature pair for volumes \ 
vol! and vol2 and get a cuirvature pair corresponding to volumes in inside a \ 
box with vl, v2 as its lower left corner*) 

(bright now this method needs to dyanically choose its changing size to avoid \ 
infinite loops*) 

(♦Returns a two element list of curvatures*) 
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(*vl, v2 should correspond to the point in the lower ri^t of the rectangle \ 
containing VolBublA[kl, k2] , VolBub2ADtl, k2] 

the curvatures kl, k2 should be a point in the middle of the rectangle*) 



NextCurvaturePairinArray [ kl_Real , k2_Real , vl_Real , v2_Real , changel_Real , \ 

change2_Real , scalingFactorl_Real , scalingFactor2_Real] \ 

: =Module [{curvatureComponentll , curvatureComponent 12 , curvatureComponent21 , \ 
curvatureComponent22, start ingVolumel , startingVoluiiie2 , volumeComponent 11 , \ 
volumeConiponentl2, volumeComponent 21 , volumeComponent 22, alpha, beta, \ 
volumeOneTarget , volumeTwoTarget , curvaturePair , volumeChangeMatrix , 

startingVoltoTargetMatrix, alphaBetaHatrix , coxmter, \ 
oneBoxSizeOver Array , oneBoxSizellpArray} , 

(♦adjustment arrays*) 
oneBoxSizeOverArray = 

MakeChangeOneBoxSizeOver [kl , k2, 
VolBublA[kl,k2] , changel, scalingFactorl] ; 

oneBoxSizeUpArray = MakeChangeOneBoxSizeUp[kl, k2, VolBub2A[kl, k2] , \ 
change2, scalingFactor2] ; 

(♦Declaration of initial curvaturevectors*) 

curvatureComponent 11 = oneBoxSizeOverArray [ [1] ] ; 

curvatureComponent 12 = k2 ; 
curvatureComponent 21 = kl; 

curvatureComponent 22 = oneBoxSizeUpArray[[l]] ; 

count er=0; 

(*Print ["cll = " , curvatureComponentll] ; 
Print ["cl2 = " , curvatureComponent 12] ; 
Print ["c21 = " , curvatureComponent 21] ; 
Print [" c22 = " , curvatureComponent 2 2] ; *) 



(^declaration of initial volumevectors*) 

startingVolumel = VolBublA[kl, k2] ; 

startingVolume2 = VolBub2A[kl, k2] ; 

volumeComponent 1 1 = VolBublA[ 

curvatureComponentll , curvatureComponent 12] ; 

volumeComponent 12 = VolBub2A [curvatureComponent 11, \ 
curvatureComponentl2] ; 

volumeComponent 21 = VolBublA[ 

curvatureComponent 21 , curvatureComponent 22] ; 

volumeComponent22 = VolBub2A [curvatureComponent 21 , \ 
curvatureComponent 22] ; 

(*Print ["vll = " , volumeComponentll] ; 

Print ["vl2 = " , volumeComponent 12] ; 

Print ["v21 = " , volumeComponent 21] ; 

Print ["v22 = ", volumeComponent22] ; *) 

(^declaration of target Volumes*) 
volumeDneTarget = vl + ,5*changel; 
volumeTwoTarget = v2 + .5*change2; 

(♦Print ["volumeOneTarget = " , volumeOneTarget] ; 
Print ["volumeTwoTarget = " , volumeTwoTarget] ; 

•) 

(♦declaration of Matricies*) 

startingVoltoTargetMatrix^ {{volumeOneTarget - 

StartingVolumel} , {volumeTwoTarget - startingVolume2}> ; 

(♦Print ["first test = 

" , IsPointinBox [volumeComponentll, volumeComponent 12 , vl+changel, \ 

v2 , change 1 , change2] ] ; 
*) 

If [IsPointinBox [volumeComponentll , volumeComponent 12, vl, v2, changel, 
change2] , 



(♦Print ["vll = " , volumeComponentll] ; 
Print ["vl2 = " , volumeComponent 12] ; 
Print ["changel = " , changel] ; 
Print ["change2 = " , change2] ; 

•) 

curvaturePair = {curvatureComponentll , curvatureComponent 12)-, 
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(♦Print ["2nd test = " , IsPointinBox [voluiiieComponent21 , \ 
voliimeComponGnt22 , vl+changel , v2 , change 1 , change2] ] ; 

Print ["v21 = " , volumeComponent21] ; 
Print ["v22 = " , volumeComponent22] ; 
Print ["changel = " , changel] ; 
Print ["change2 = " , cha[ige2] ; 

•) 

If [IsPointinBox [voluiiieComponent21 , voliimeComponent22 , vl , v2 , changel , 
change2] , 

(*Print["case II"] ; 

Print ["v21 = " , voIumeComponent21] ; 
Print ["v22 = " , voIumeComponent22] ; 
Print ["changel = " , changel] ; 
Print ["change2 = ", change2] ;*) 

curvaturePair = ■[curvatureComponent21 , curvatureComponent22>, 

(*Print ["while loop test", \ 
! IsPointinBox [volumeComponent 11 .volumeComponent 12, vl+changel , v2, changel , 
change2] ] ; 

*) 

voIumeChangeMatrix = -t-[voIumeComponentll - start ingVolumel, \ 
volumeComponent 12 - start ingVoIuin82}, {volumeComponent 21 - 

St art ingVolumel , volumeComponent22 - start ingVolume23-} ; 

(♦Print ["voIumeChangeMatrix = " , voIumeChangeMatrix/ /MatrixForm] ; 

•) 

alphaBetaMatrix = Inverse [voIumeChangeMatrix] . 

startingVoItoTargetMatrix; 
(♦Print ["alphaBetaMatrix = " , alphaBetaMatrix/ /MatrixForm] ; 

•) 

alpha = alphaBetaMatrix [ [1, 1] ] ; 
beta = alphaBetaMatrix[[2,l]] ; 

curvatureComponentll = CcurvatureComponentll-kD^alpha + \ 
(curvatureComponent21-kl)+beta +kl ; 

curvatur eComponent 12= (curvatureComponentl2-k2)^ 
alpha + (curvatureComponent22-k2)^beta + k2; 

(♦Print ["inside while loop"]; 

Print ["cll = " , curvatureComponentll] ; 
Print ["cl2 = ", curvatureComponent 12] ; 
Print ["alpha = ", alpha]; 
Print ["beta = ", beta];*) 

volumeComponent 1 1 = VoIBublA[ 

curvatureComponentll , curvatureComponentl2] ; 
volumeComponent 12 = VolBub2A[ 
curvatureComponent 1 1 , curvatureComponent 12] ; 

(*Print ["vll = " , voIumeComponentll] ; 
Print C"vl2 = " , voIumeComponentl2] ;♦) 

]; 

curvaturePair = {curvatureComponentll , curvatxireComponentl23-; 
(♦Print ["after loop"] ; 

Print ["vll = " , voIumeComponentll] ; 

Print ["vl2 = " , volumeComponent 12] ; 

Print ["changel = " , changel] ; 

Print ["change2 = " , change2] ; 

Print ["2nd 

test = " , IsPointinBox [voIumeComponentll , \ 
volumeComponent 12 , vl+changel , v2 , changel , change2] ] ; *) 



]; 
(♦If [! 

IsPointinBox [VolBubl [curvatureComponentll , curvatureComponent 12] , \ 
VolBub2 [curvatureComponentll , curvatureComponent 12] , vl+changel, v2, \ 
changel, change2] , 

curvaturePair = {-1, -1}; 

]; 
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curvaturePair 

]; 



C*This retums a lower bound on the positive part of the Hutchings function \ 
for a given volume pair*) 

LowerBoundOnPosHutchingFunct ion [voll_Real , vol2_Real , errorl_Real , \ 
error2_Real] := 

2*AreaSphere [Curvaturef romVolSingle [voll/2 , errorl] ] + 
AreaSphere [Curvaturef romVolSingle [vol2 , error2] ] + 
AreaSphere [Curvaturef romVolSingle [voll+vol2, error2] ] ; 



C*This returns an array of areas is indexed by its volume pair *) 
ArrayBuilderForSingAreasandCurvatures [vMin_Real , vMax_Real , boxSize_Real , \ 
error _Real] :=Module[ 

{AreaArray , volumeArray , n , g} , 

g [i_Integer] = AreaSphere [Curvaturef romVolSingle [vMin+ i*boxSize , \ 
error]] ; 

n[i_Integer] =VolSphere [Curvaturef romVolSingle [vMin+ i*boxSize, error]] ; 

AreaArray= Array [g, Ceiling [CvMax-vMin)/boxSize] , 0] ; 
volumeArray=Array [n, Ceiling [CvMax-vMin)/boxSize] , 0] ; 

{volumeArray , AreaArray} 

]; 



(♦This retums an array of areas is indexed by its volume pair *) 
ArrayBuilderForSingAreasandVolumes [vMin_Real , vMax_Real , boxSize_Real , \ 
error _Real] :=Module[ 

{AreaArray , volumeArray , n , g} , 

g[i_Integer] = AreaSphere [Curvaturef romVolSingle [vMin+i*boxSize , 
error]] ; 

n [i_Integer] =vMin+i*boxSize ; 

AreaArray^ ArrayLg, Ceiling[CvMax-vMin)/boxSize] , 0] ; 
volumeArray=Array [n , Ceiling [(vMax-vMin)/boxSize] , 0] ; 

{volumeArray , AreaArray} 

]; 



(*This returns an array of areas is indexed by its volume pair *) 
ArrayBuilderForSingAreas [vMin_Real , vMax_Real , \ 
boxSize_Real , error _Real] : =Module [ 
{AreaArray, g}, 

g [i_Integer] = AreaSphere [Curvaturef romVolSingle [vMin+i*boxSize , \ 
error] ] ; 

AreaArray= Array[g, Ceiling [CvMax-vMin)/boxSize] , 0] ; 
AreaArray 

]; 



(*This code fills an array with curvature pairs. It moves on to the next \ 
region after verifying that 2 the area of the double bubble enclosing volumes \ 
v,w is less thein the concave part of the Hutchings function for v,w*) 
ArrayFillingPrcof [vMin_Real , wMin_Real , wMax_Real , rectangleHeight_Real , \ 
rectangleWidth_Real , startingCurvaturel_Real , startingCurvature2_Real , \ 
adjustmentMainVol_Real , adjustment SecondVol_Real] : =Module [{smallerVolume , \ 
largerVolume,nextRowStartingCurvatures, curvaturePair, f ailSaf e, 

nextRowStartingPosit ion, counter, insideCounter , f ailingVolumeV, \ 
singAr easArray , singAreasArrayStart}, 

failSaf e=True; 
count er=0; 

insideCounter =1 ; 

singAreasArrayStart= (vMin-rectangleWidth) /2 ; 

singAreasArray= ArrayBuilderForSingAreas [(vMin-rectangleWidth)/2, \ 
2*CwMax+rectangleWidth) , rectangleWidth/2, rectangleWidth/2] ; 
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Print ["Array complGted, evaluating Hutchings Function"] ; 
(*Print ["singAre as Array =" , singAreas Array] ; *) 

(♦Print ["hi"] ; 
Print [f ailSaf e] ; *) 

C*nextRowStartingPosition keeps tract of the x position that \ 
starts a new row 

initially it is set to and then set to correct position*) 
nextRowStartingPosition = 0.0; 
nextRowStartingCurvatures = curvaturePair ; 



curvaturePair = {startingCurvaturel , startingCurvatureS}; 



smallerVolume = vMin; 
largerVolume=wMin+rectaiigleHeight ; 

(♦Print ["first test", largerVolume wMax +rectangleHeight ftfe \ 
failSafe] ;*) 

While [largerVolume wMax +rectangleHeight failSafe, 
(♦While the vl is less than v2*) 

VHiile [smallerVolume largerVolume +rectangleWidth failSafe , 

failSafe = IsPointinBox [VolBubV [curvaturePair [ [ 
1] ] , curvaturePair [ [2] ] ] , 

VolBubW[curvaturePair[[l]] , curvaturePair [[2]]] , smallerVolume, \ 
largerVolume , 2^rectangleWidth , 2^rectangleHeight] ; 



If [IfailSafe, 



Print ["( 

VolBubV [curvaturePair [ [1] ] , curvaturePair [ [2] ] ] , " , " , 
VolBubW [curvaturePair [[1]] , 

curvaturePair [ [2] ]] , ") is not in the 

box defined by (", smallerVolume,", ", largerVolume,") \ 
and (" , rectangleWidth, " , " , rectangleHeight , ") " ] ; 



] ; 

(♦marker*) 

(♦this sets up moving one row up^) 

If [ nextRowStartingPosition \ [Equal] && .85 * largerVolume \ 
smallerVolume+rectangleWidth, 

nextRowStartingPosition = smallerVolume; 
nextRowStartingCurvatures = curvaturePair ; 

]; 

(♦Print["hi! !!"];♦) 
(♦If [insideCounter <= 5, 

Print ["target volume pair = smallerVolume, ",", largerVolume, 
">"] ; 

Print ["actual volume 

pair = \ 

VolBubV [curvaturePair [ [1] ] , curvaturePair [ [2] ] ] , " , " , \ 
VolBubW [curvaturePair [[1]] , curvaturePair [ [2] ] ] , "}"] ; 

Print ["curvatures are " , curvaturePair] ; 

insideCounter++ ; 

Print ["failSafe ", failSafe]; 

];*) 

(♦if there is a problem with the curvature function or the \ 
hutchings function is negative then the program should fail^) 

(♦Print ["v/2 from array =" , \ 
singAreasArray [ [Floor [ ( (smallerVolume-rectangleWidth) /2-singAreasArrayStart ) / (\ 
rectangleWidth/2)+l] ] ] ] ; 

Print ["index =" , \ 

Floor [( (smallerVolume-rectangleWidth) /2-singAreasArraySt art)/ (rectangleWidth/\ 
2)+l]]; 



Print ["v/2 not from array = " , 

AreaSphere [Curvaturef romVolSingle [( 
smallerVolume-rectangleWidth) /2 , \ 

rectangleWidth/2] ] ] ; 



50 



J. CORNBLI BT AL. 



Print ["w from array =", 
singAreasArray [ [Floor [ 

2+(largerVolume-rectangleWidth-singAr©asArrayStart) / C\ 
rectaiigleWidth/2)]]]] ; 



Print ["index =" , \ 

Floor [2+ (larger Volume-rectangleWidth-singAreasArraySt art ) / (rectangleWidth/2) ] ] \ 



Print ["w not from array = " \ 
, AreaSphere [CurvaturefromVolS ingle [larger Volume- 

rectangleWidth, rectangleWidth/2]]] ; 

Print ["v+w from array =", \ 
singAreasArray [ [Floor [1+ (smaller Volume+largerVolume-rectangleWidth- 
singAreasArrayStart) / (rectangleWidth/2)]] ] ] ; 

Print [" index =" , \ 

Floor [1+ (-rectangleWidth+smallerVolume+largerVolume-singAreasArrayStart )/ (\ 
rectangleWidth/2)]] ; 

Print [ 

"v+w not from array = " \ 
, AreaSphere [CurvaturefromVolS ingle [smallerVolume+largerVolume-rectangleWidth, \ 
rectangleWidth/2] ] ] ; 

Print ["v + w =" , smallerVolume+largerVolume-rectangleWidth] ; 

Print ["v = " , smallerVolume] ; 

Print ["w = larger Volume-rectangleWidth] ; 

*) 

If [curvaturePair[[l]] < 1 \ 
I 1 2*singAreasArray [ [Floor [( (smallerVolume-rectangleWidth) /2-\ 
singAreasArrayStart)/(rectangleWidth/2) +1] ] ] +singAreas Array [ [ 

Floor [2+ (largerVolume-rectangleWidth-\ 
singAreasArrayStart)/(rectangleWidth/ 

2)] ] ] +singAreasArray [ [Floor [l+(smallerVolume+\ 
largerVolume-rectangleWidth-singAreasArrayStart) /(rectangleWidth/2) ] ] ] + - \ 
2*AreaDblBubble [curvaturePair [ [1] ] , curvaturePair [ [2] ] ] < , 

Print ["curvaturePair =" , curvaturePair] ; 

Print ["hut chings function = \ 
2* SingAreasArray [ [Floor [ ( (smallerVolume-rectangleWidth) /2- 

singAreasArrayStart)/ (rectangleWidth/2)+l] ] ] +singAreas Array [ [\ 
Floor [2+(largerVolume-rectangleWidth-singAreasArrayStart)/ (rectangleWidth/2)] ] \ 
] +singAreasArray [ [Floor [1+ (smallerVolume+largerVolume-rectangleWidth-\ 
singAreasArrayStart) / (rectangleWidth/2) ]] ] - 2*AreaDblBubble [ 
curvaturePair [ [1] ] , curvaturePair [ [2] ] ] ] ; 

failSafe = False; 

f ailVolumeV =smallerVolume ; 



Print ["upper 

bound on dbl bubble", 2*AreaDblBubble [curvaturePair [ [1] ] , \ 
curvaturePair [ [2] ] ] ] ; 

Print["lower bound on concave" .singAreasArray [ [ 

Floor [ ( (smallerVolume-rectangleWidth) /2-\ 
SingAreasArrayStart) / CrectangleWidth/2) +1] ] ] +singAre as Array [ [Floor [2+(\ 
largerVolume- 

rectangleWidth-singAreasArrayStart)/ CrectangleWidth/2)] ] ]+\ 
SingAreasArray [ [Floor [l+( smallerVolume+largerVclume-rectangleWidth- 
singAreasArrayStart)/(rectangleWidth/2)]]]] ; 

]; 

(♦Print ["inside while loop"];*) 

curvaturePair = NextCurvaturePairinArray [curvaturePair [[1]] , \ 
curvaturePair [ [2] ] , smallerVolume+rectangleWidth , largerVolume , 
rectangleWidth, rectangleHeight , adjustmentMainVol, \ 

adjustmentSecondVol] ; 

(*Print ["upper bound on dbl bubble", 2*AreaDblBubble \ 
[curvaturePair [ [1] ] , curvaturePair [ [2] ] ] ] ; 

Print ["lower bound on \ 
concave" ,LowerBoundOnPosHutchingFunct ion [smallerVolume , largerVolume, \ 
rectangleWidth, rectangleHeight]] ; 
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smallerVolume = smallerVolume+ rectangleWidth; 



] ; 



If [failSafe, 

(♦Print ["smallerVolume" , smallerVolume] ; *) 
smallerVolume = nextRowStartingPosition; 

nextRowStartingPosition =0.0; 
count er++; 
insideCounter= 1 ; 



If [counter \ [Equal] 100, 

Print ["w at ", largerVolume] ; 
counter=0 ; 

]; 



(♦Print ["SETTING the curvature pair"]; 

Print ["Volumes before {" , \ 
VolBubV [nextRowStartingCurvatures [ [1] ] ,nextRowStartingCurvatures [ [2] ] ] , " , \ 
" .VolBubW [nextRowStartingCurvatures [ [1] ] , 

nextRowStartingCurvatures[[2]]] , "}" ]; 

Print ["target box { (" , smallerVolume " , " , \ 
smallerVolume+rectangleHeight , ") X C 

" , (largerVolurae+rectangleHeight) , " , " , \ 
ClargerVolume+2*rectEingleHeight) " ) >"] ; 

*) 

largerVolume = largerVoluine+ rectangleHeight ; 

While [ ! IsPoint inBox [VolBubV [curvaturePair [ [1] ] , 

curvaturePair [ [2] ] ] , VolBubW [curvaturePair [ [1] ] , curvaturePair [ [ 
2] ] ] , smallerVolume , largerVolume , 2+rectangleWidth , \ 
2^rectangleHeight] , 

curvaturePair = NextCuirvaturePairinArray [curvaturePair [ [1] ] , \ 
curvaturePair [ [2] ] , smallerVolume , largerVolume+rectangleHeight , \ 
rectangleWidth, rectajigleHeight , adjustmentMainVol , adjustmentSecondVol] ; 
] ; 

(♦Print ["VOLl =" ,VolBubV[curvaturePair[[l]] , curvaturePair [ [2] ] ] ] ; 
Print["V0L2 =",VolBubW[curvaturePair[[l]] , curvaturePair [ [2] ] ] ] ; 
♦) 



failSa±e=IsPoint inBox [VolBubV [curvaturePair [ [1] ] , curvaturePair [ [2] ] ] \ 
, VolBubW [curvaturePair [ [1] ] , curvaturePair [ [2] ] ] , smallerVolume , largerVolume , \ 
2^rectangleWidth, 2^rectangleHeight] &ftsmaller Volume/ largerVolume . 85 ; 

If [IfailSafe, 

Print ["error at end ( ", VolBubV [ 

curvaturePair [ [1] ] , curvaturePair [ [2] ] ] , " , 
" , VolBubW [curvaturePair [ [1] ] , curvaturePair [ [2] ] ] , " ) 
is not in the box defined 
by ( " , smallerVolume , " , " , largerVolume , " ) 
and (" , rectangleWidth, " , " , rectangleHeight , ") " ] ; 



(♦Print ["test " , IsPointinBox [VolBubV [curvaturePair [ [1] ] , 
curvaturePair [ [ 

2]]] , VolBubW [curvaturePair [[1]] , curvaturePair [ [2] ] ] , smallerVolume, 
largerVolume , rectangleWidth , rectangleHeight] ] ; ♦) 

]; 

]; 

(♦Print ["while loop test ", smallerVolume largerVolume && failSafe, \ 
failSafe] ; 

Print ["first part of test", smallerVolume largerVolume] ; 



If [failSafe, 
Print ["done"] , 
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Print ["failed at (" , smallerVolume , " , " , largerVolume , ") "] ; 

Print ["ratio =" , smallerVolume/largerVolume] ; 

]; 

Print ["failsafe = " , f ailSaf e] ; 
]; 
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